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PREPACE 


The Department of Aeronautics, Naval Postgraduate School 
offers several courses at different levels on the dynamics of 
flight vehicles. These courses include studies of fixed wing 
Stability, rotary wing stability, and the stability of 
missiles. In addition, there are general courses that cover 
the dynamics of bodies in six degrees of freedom, including 
the development of the equations of motion. 

There is currently no signle resource available that 
covers the material for the airplane and missile courses, as 
well as the introductory dynamics classes, and, as a result, 
the instructors must make use of large numbers of locally 
prepared handouts and/or assign several textbooks, only 
portions of which will be used. 

In order to provide a solution to this problem, the 
following material has been developed for use in the dynamics 
classes in the Department of Aeronautics. Although this 
material is, in places, of greater depth than what would be 
covered in the courses, it is felt that the value of this 


material is enhanced by its degree of thoroughness. 
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CHAPTER l 


DYNAMICS OF MOTION 


1.1 INTRODUCTION 

This chapter is intended to be a formal review of 
mathematical and physical concepts that will be used to 
develop the theory of dynamic stability of ~flight vehicles. 

The starting point is the basic laws of Newtonian 
Mechanics, with a useful application to spring-mass-damper 
systems. We shall find later how the motion of any 
conventional flight vehicle is similar to those systems. 

We then proceed to rigid body equations of motion, 


described in non-inertial reference frames. 


1.2 NEWTON'S LAWS OF MOTION 

In any system of units, basic quantities and units are 
related through Newton's laws of motion. 

NEWTON'S 1ST LAW OR THE LAW OF COORDINATE SYSTEMS: 

If there is no force acting on a particle, its velocity 
remains constant. 

This law is often called the law of inertia and applies 
to inertial frames of reference. 

NEWTON'S 2ND LAW OR THE LAW OF DYNAMICS: 

The time rate of change of the linear momentum of a 


Particle, is proportional to the net force acting on it. 


14 





The vector quantity 
pet) = m(t)- v(t) (ai) 
is defined as the linear momentum of the particle, where m(t) 
is the inertial mass of the particle and v(t) is its velocity 
relative to the reference frame in question. 


The 2nd law is mathematically expressed as 


c t dept) 
nee oa 
9 © (1-2) 
Or 
= 
Fined ae. ae = Panty. D | 
Cl) 
If the mass m remains constant, we can write: 
a { d = 
eee — t 
F med a. oe, v (t) (1-4) 
Or 
-—p { —_ 
=—_ ma (t 
Feet gma (1-5) 


where a(t) is the acceleration caused by the net applied 
force eae 

The proportionality factor is expressed as 4 and its 
value depends on the system of units being used. 

In the English Engineering system, one pound force 
(l-lbf£) will cause to one pound mass (l-lbm), an acceleration 
of 32.174 ft/sect Thus substituting into Equation 1-5 

L-Ihf = - (Ibm) 32.174 (4/ sec?) a 


10S 





Hence 
9. = 32.174 Ibm ft /\bf -sec? (1-7) 

A few other units used for velocity, are given below with 

their corresponding conversion factors: 
1 mile per hour (lmph) = 1.4666 ft/sec 
1 knot = 1 nautical mile per hour = 1.6878 ft/sec 

NEWTON'S 3RD LAW OR THE LAW OF MUTUAL FORCES: 

If one body exerts a force on a second body, the latter 
exerts an equal but opposite force on the first. 

This law allows us to define the gravity force, or the 
weight of a body, that has a given mass (m). It is the force 
that results due to the earth's attraction. 

Besides the three fundamental laws of mechanics, the 
Newton's Universal Gravitational Law states that two inertial 
masses ml and m2, Separated by a distance R are subject to an 


attraction force according to the equation: 


am, MM 2 


F=QeG a: (1-8) 





where G, is the gravitational constant, which in the EE 


system is G = 6.32«10. !bf-ft? /ibm? 


5 
If now m, = Mp = 1.32 X 10°° lbm (mass of the earth), 
then near earth's surface where R = Re = 2.99 X Tale ft 


(radius of the earth), we can write: 


G Me 
WN = cle sya tM (1-9) 
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or 

Vance (1-19) 
where g = g. GMe /Re = 32.174 ft/sec*, called acceleration of 
gma. ty. 

By Equation 1-18 we can see that l-lbm, subject to the 
earth's gravitational field, experiences a force l-lbf and 
hence it accelerates with 32.174 ft/sec’, near earth's 
Surface. 

We can also define another unit of mass, the slug, as the 
mass that experiences an acceleration l-ft/sec*, if subject 
mo a force of 1-lbf. 

By this definition 

1 slug = 32.174 lbm (1-11) 

Consequently, we can derive the density (¢), as the mass 


per unit volume usually measured in slugs/ft. 


muss (™m ) 
“yolume CV) _ (1-12) 


density (p) = 
People in the field use the (gq), as a short hand notation 
for force or acceleration. A 5-g force simply means a force 


5 times the body's weight. 


1.3 ANGULAR MOMENTUM 
A special significance to the 2nd law is given by the 
angular momentum approach derived simply by crossing the 


radius vector r ; on to the 2nd law. 


If 





5 ag = x Pek 
dt Cie) 
We define 
Gmet = rx co (1-14) 


the net torque or moment of the force about the origin and 


— 


h=rxp (1-15) 
the angular momentum about the origin. 


Differentiating Equation 1-15 with respect to time, we 


get 


cit dt alt (1-16) 


The first term is zero, because 





dr — = ~ - 

——— = — “yn 4 

a xD ux(mv) (uxv) =O een 

Figure l-l 
Moment of Force About the Origin 
So we finally get 

dh _ 2 1-18 
at = Conct ( ) 
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which tells us that the time rate of change of the angular 
momentum of a particle, equals to the net torque or moment 
acting on it. 

If no torque is acting on a particle, its angular 
momentum 1s conserved, likewise its linear momentum is 
conserved, if no force is acting on it. 

A simple example 

For a pendulum that swings in a vertical plane, having a 
mass (m), the acting forces are, its weight W = mg and the 
tension of the rod T. The rod is assumed to have negligible 


mass. 





Figure 1-2 


A Simple Pendulum 


If 9 is the angle between the rod and the vertical axis, 
s, is the arc length corresponding to the angle 9 and l, is 
the length of the rod, a direct application to the 2nd law 


(Equation 1-5), gives: 


eo 





-masin@ = ms | (1-19) 
and equilibrium condition in the direction of the rod, gives: 
-mag cos © = if (1-29) 
Equation 1-19, can be written 
S$ +gsin@ =O (1-21) 
If 9 is small, we eon as sume 
0 = a = 5 = 10 eae 
and 
6 = sino ka) 
With this assumption, Equation 1-21, can be written 
6 + + @ =0 
(1-24) 
Equation 1-24, is a linear differential equation, with 
constant coefficients, which describes the motion of the 
pendulum. 
The same result can be derived, if we use the angular 
momentum approach. 
The angular momentum of the pendulum, about the pivot, is 
Simply ({), cross the linear momentum (ms). The time rate of 
change of the angular momentum, has to be equal to the net 


external torque. 


d ; 
Bee: —— © 
2 (2x ms ) £ ang sin oe 
Since 4 is constant (dl/dt = 8), we can write 
Qans = ~ €mgsinO (1-26 ) 
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Or 
6+ ma cine =. (1-27) 
or 


le q 
=o = 
0 +5 e (1-28) 
The minus sign in Equations 1-19 and 1-25, indicates that 
if the angle 9 decreases, the acceleration $ in Equation 1-19 


and the angular momentum in Equation 1-25, increase. 


1.4 UNDAMPED HARMONIC OSCILLATOR 

In this section and in the following section, we will 
derive the equations and the solutions to the spring mass 
damper systems. 

The undamped oscillator is simply a spring mass system 
without friction, shown in Figure 1-3. 

Let k, denote the spring stiffness, m, the mass of the 
body, x(t), the displacement of the body from the equilibrium 
position. Notice that the weight of the body stretches the 
Spring permanently, so that xe is the equilibrium position, 


after the body has hanged to the spring. 


Xe k = 
ay hs 
wm 

Prgunre:.—s 


Spring Mass System 
21, 





The restoring force according to Hook's law, is 
=x (1-29) 
Then according to the 2nd law 
mx =-Kx (1-39) 
Or 
(1-31) 


As an example we will take 


m= 18 lbm (1-32) 
k = 40.5 1bf/ft (1-33) 
If we let 
2 k 
Wy = Ye Pe (1-34) 


where Wy is the so called natural frequency of oscillations, 
Equation 1-31 becomes 
X + Wy Xx =O (1-35) 
This is a linear differential equation with constant 
coefficients. Its characteristic: equation is clearly 
q+, = 0 (1-36) 
which has characteristic roots 
Q = 2 hy (1-37) 
The most general solution of Equation 1-35 can be 


expressed as 


(1-38) 


a2 





where A, and A_ are constants depending on the initial 


conditions. 


Expanding the exponential terms of Equation 1-38, we 


obtain 
x(t) = o& coswyt + B SIN Wnt (1-39) 
where 
x= A,+A_ (1-40) 
and 
Oi tiene) (1-41) 


If we now express 
sin Wat = cos (Wnt + 90°) (1-42) 
we can give the following useful form to our solution 
x(t) = Xmay cos (Wnt + Oe) (1-43) 
where Xmqgx is the maximum deviation from the equilibrium 
position and %&, the initial phase. 
We can relate the quantitieS Xm, and % , with the 
initial conditions, as follows: 
At time t=@8, the initial displacement x, , is given by 
Equation 1-43, 
4s Morne: CoS ee (1-44) 
aiamtne inttial velocity X, , is given by differentiating 
Equation 1-43 with respect to time and evaluating at t=6. 
Se Saye De Gin lee (1-45) 
As an example we will take 
X. = 4 ft (1-46) 


Xo 


32 ft/sec (1-47) 
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From Equations 1-44 and 1-45, the initial phase and the 


maximum displacement can be found to be 





a 

9, =+an (=) (1-48) 
sae, 
2 oe 

Xmax = ( Xe + we (1=49) 


Note: In Equation 1-48, keep the sign, so that the 
quantity in brackets will give the right quadrant for the 
angle 0, e 


For our example 


K 40.5 
Wy = g. on 32.174 » ie 79.39 sec” 








(1-56) 
7 =| S27 : y+ 
©, = tar’ (—"——) = 436.76" = 2.59 rad ied 
Monae = H(d6e 1024 )" oa 
a 72.39) ~~ ~ (1-52) 
so that the response is 
x(1) = -5.49 cos( 8.54 4 + 2.39 pad) (1-53) 


1.5 DAMPED HARMONIC OSCILLATOR 
Consider the spring mass damper system shown in 
Figure 1-4. Assume that the friction force exerted by the 


damper, is a linear function of the mass velocity. 
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kt ee a 
<A) 4 Fs 
T 7 


Figure 1-4 


Spring Mass Damper System 


Let k the spring stiffness,-m the mass of the body, x(t) 
the displacement of the body from the equilibrium position 
Xe, aS it was defined in section (1.4), and c, the 
coefficient of friction. 

The restoring force is again 

Pa @ =18% (1-54) 


and the friction force is 


Ine (1-55) 


Then according to the 2nd law 


am X Se Ge Xx (1-56) 


Or 
(1-57) 


Z5 





As an example we will take 


m= 18 lbm (1-58) 
k = 40.5 lbf/ft (1-59) 
If we let 
2 k 
Wn = Y. = (1-60) 
and 
T= ae 
mon: adc a (1-61) 


where 5 is the so called damping coefficient, then the 
Equation 1-57 is written 
% + 25Wyx + Wn X =O (1-62) 
This is a linear differential equation with constant 
coefficients. Its characteristic equation is clearly 
9? + 23w_ 1 + We = 0 (1-63) 
Three cases arise in obtaining solutions to the 
differential Equation 1-62. 


AS an example, we will take 





Xe = 4 ft (1-64) 
% = 32 ft/sec (1-65) 
: sme Beli 
As an example we will take 
c = 952 Ibf /tt- sec (1-66) 
so that : 
aloe | (1-67) 


In this case, the characteristic equation has a double 


root 


aaa, (1-68) 
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and the most general solution is 
x(t) = (A, + At] oo unt (1-69) 
where A, and A_ are constants determined by the initial 
conditions as follows 
A, = Xo (1-76) 
and 
A_ = Xo + XeWn Cleat 
so that the solution in terms of the initial conditions can 


be written 
~wnt 


x(t) = [Xo (dt+unt) + xX +] © (i=7y 
For our example 
x ({) =-[4 (1+ 8514) + 324] :o ed Si 
Or 
x(t) = (fet, Gee Can (1-74) 
Case 2: Overdamped Motion (5>1) 
As an example we will take 
c = 14.28 lbf/ft-sec™ (1-75) 
so that : 
5 = 9. ee (1-76) 


In this case the roots of the characteristic equation are 


real and negative 
)" 


4 = -5w, + Wy (5°-1 (T=77) 
and the most general solution is 
~5wnt ; Wal 5*-4 t -Wn¥ 2-4 t | 
i= C | A, @ +A (T=78) 
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where A, and A_ are determined by the initial conditions as 
follows 


ees NaN rue7s) 


% = = Bw (ALtA.) twy V541 (A,-A_) (1-88) 
which yields 


((FT+5) x + 
| pe ene SS ai 


9 (S41 (diets 

, , LRT -8)4~ 
- 25-5 (1-82) 

For our example 
A, = 6.36 ft (1-83) 
A_ = -2.36 ft (1-84) 
and 

x(1} = go (6.36 oa OESe eae : } (1-85) 

or / 
Mes etc e Sheee *** (1-86) 


Case 3; Underdamped Motion (5<]) 


As an example we will take 


c = 4.76 lbf/ft-sec"™ (1-87) 
so that 
Cc 
One OS 
5 = 9. Den irs (1-88) 


In this case the roots of the characteristic equation are 


imaginary with negative real parts. 


g = -5Wn tj Wy 11-8 (1-89) 
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If we define 
Wy = Wn (4-5? (for 5<1) (1-96) 
the so called damped frequency of oscillations, 
2 -S0n tj wy (1-91) 


and the most general solution is 


Sn t jw, t -jW,t 
xt) =e (Ae 7 +A ) (1-92) 
Expanding the exponential terms we obtain 
-Swnt 
xt) =e — (acosayt +8 sin w,t) (1-93) 
where 
= A, + A- (oa 
and 
B= 5 (A, -A-) (1-95) 
The solution can be formulated as 
Sunt 
Vitex c cos (ut +O.) (1-96) 


where X and 9 are constants determined by the initial 
conditions as follows: From Equation 1-69 we obtain 
nu 
= x(t) = X cos Cwyt+ 9) (1-97) 


Differentiating Equation 1-97 with respect to time we obtain 


— | Sw, x(t) + x(t) | =-Xw, sin (w, +t +0.) (1-98) 


For t=@, Equations 1-97 and 1-98 become respectively 


= X cos O,5 (1-99) 
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and 


SW Xo t %y = -XwWy sin O, (1-198) 


So the constants are determined as 





a Ge 
9, = tan | Sun Xo + nee | (1-101) 


and 


xe + “ (SWn Xo & %) |" 
; | ae wa? (1-192) 
Notes In Equation 1-191, keep the sign, so that the 


quantity in brackets will give the right quadrant for the 


angle oa 
For our example 
/ : 
wy = 851 (1-057)? = 7.37 sec! (1-103) 
Semen 02) = 25lherad (1-104) 
X = + 7.76 ft (1-195) 
so that the response is 
~4. 26+ 
x(t) =-%76 e cos (F7.3%+ + ZI rad) (1-186) 


The roots of the characteristic equation can be plotted 
in the complex plane as shown in Figure 1-5. 

In cases 1 and 2, the damping is so large that no 
oscillations take place and the mass simply returns gradually 
to the equilibrium position x=@. In those cases the roots 
are real and negative. 

In case 3, the damping has been reduced, so that 
oscillations about the equilibrium position, do take place, 


the amplitude of which dies out as t7®, 
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Wa © Wy (1-59) 





* 


~~ ee = & ~=~ = @ ae a= ao 


%- 


‘Figure 1-5 


Plot of the Characteristic Roots in the Complex Plane 


Figure 1-6 shows the responses in the above three 
described cases where it is assumed that the initial 


displacement is x. , and the initial velocity is %=0. 





critically damped motion 51 

overdamped motion: 5>1 
underdarnped motion 5<i 
envelope 


head 
= 
=~ « 


Figure 1-6 


Response of Displacement (x) of a Spring Mass Damper System 


Sal 





Notice that in the case where the damping 5 is zero, we 
have a spring mass system, without friction, described in 
section 1.4. In this case, the roots are imaginary and w, is 
the same as Wn (Figure 1-5). 

It is clear that the system's physical constants (in our 
case 5,m,c,k), uniquely define the roots of the 
characteristic equation. Those roots are given the name: 
normal modes because one can determine from those roots, the 
corresponding characteristic motion (in our case, undamped 
Critically damped, overdamped, underdamped). 

We have seen so far how an harmonic damped oscillator 
responds, if initially disturbed from its equillibrium 
position. Then Knowing the initial conditions we can always 
find the response of the displacement or velocity versus 
time. 

Now we are interested in what happens if the oscillator 
is driven by an external input force F(t). 

In this case the solution will consist of two parts: 

x(t) = Xpom 4) + Xpara () (1-197) 
where Xhom(t) is the homogenous part of the solution 
discussed so far, i.e. considering the input force zero and 
puri (4) is the particular solution due to this particular 
applied force, i.e. considering the initial conditions zero. 

To obtain particular solutions one can use the method of 
undetermined coefficients, the method of variation of 


Parameters or the Laplace Transform method. 


a2 








In illustrating the above methods we are considering the 
following cases: 
(a) Sinusoidal Inputs 
The spring mass damper is driven bythe sinusoidal force 
F(t) =F, cosw t (1-188) 
As an example, we will take 
EGE =a COs sist t (1-189) 
The method of undetermined coefficients assumes a 
Superposition of Sinusoidal solutions of the form 
Xpart (t) = Ay coswt + AL sinwt (1-118) 
or equivalently 
X pare (4) = X cos (wt+ ) (1-111) 
where X and 9 are the coeffieients to be evaluated. 


Substituting into the differential equation 


x+25w,x +wW,x = F coswt (1-112) 
we obtain 
(wi-w*) cos(wt+Oo) - 25 wnw sin (wh +90) = 2 cos wt (1-113) 


Expanding cosine and sine and collecting terms, we get 


L (wy -w?) cos 0, - 25Wnwsin 0, | cos wt 


-| (wy -w?) sin QO, + 25WyW cos Q, | sinwt = = coswt (1-114) 
Equating coefficients 
(wr -w*) cos 8 - 23wnwsind, = -- (1-115) 
and 
( wy -w*) sin ®, + 23wyw cos @ =0 (1-116) 
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From Equation 1-116 we obtain 


Jhees a) -25 WW 


eageee (1=197) 
Multiplying Equation 1-115 by cos&®, Equation 1-116 by 


sin%, and by adding, we obtain 


Fo 


oy Sb) = cos 9, 


DS (1-118) 





The cos & is expressed in terms of tan % as 





so = _ *Fan®. Ww ae 
Adio - 
[i+ (Find yy" [ Cung-uot)* + (250m)? 
so that . 
x (1-129) 
L (we -w?)? +(23wnw)? |? 
and 
9 Sy -- -25WnW 
° oo Gene (1-121) 


For our example with the three different values of 5 


considered, the responses are: 


for 5=4 Xpare (4) = 0.28 cos (35+ +55 rad) (1-122) 
for 5=1.5  — Xgard (4) = 0.22 cos (354+5231 rad) (le12e) 
for 5=0.5 Xpard 1) = 0.36 cos (35+ + 5.82 rad) (1-124) 


The total responses correspondingly are: 


-35\t 
x(t} = (4+66.04t) e (1-125) 
+ 0.28 cos (3.54 +5.5 rad) 


~3.26t = 
x(t) = 6.36 e = 3G geo 


+ 0.22 cos (3.54 +534 rad ) 


(1-126) 
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PAG 
v= =arore pag te tierad ) 
+ 036 cos (35+ +5.82 rad ) (1-127) 


(b) Singularity Inputs 

The Laplace Transform method can be very easily applied, 
Since we are dealing with linear differential equations 
having constant coefficients. 

Appendix A gives a short Laplace Transform Table and some 
Laplace Transform Theorems that help to transform 
differential equations to algebraic. 

In illustrating the above method, we are considering unit 
impulse and unit step input forces. 

(1) Unit Impulse Response 

The spring mass damper is driven by the singularity 
input 
FE} = 0G) (1-128) 
The differential equation becomes 
% + 25Wn x + Wn x = St) (1-129) 
Neglecting initial conditions and Laplace 
transforming, we obtain: 
Xoare (S) = (st4 250ns sue) = 1 (1-138) 
Or 
X part (s) = 4 / 57+ 25 nS +Wn (1-131) 


For our example with 5=1 


X part (S) = 4 /054854)? 


-8.5it 
* par4 (4) = t @ (1-133) 


(l= 32Z) 


a5 





and the total response is 


~85i+t ~B5it 
x(t) = (4+66.04t) e@ Glee (1-134) 


with 5=1.5 


0.0525 , 0,0525 


= 4 : 72.4 = 4 G 2) ——_= 
Xart (s) = 4/574 25.535 +72.42 = 1/(5+3.26)(s+ 22.28) oe eae 





(1-135) 
-22.28+ 
Xeart (4) = 0.0525 oo ae ~ 0.0525 ae : (1-136) 
and the total response is 
~%.26t E27 284 -3.76t -22.28+t 
MIENGSS € | =.2:36 ¢ +0.0525@° 0.0525 é” (1-137) 
With 5=0.5 
Xourt (8) = 1/524 851543242 = Voi 4 26 +j £37)C 54426 -; 137) 
j 90° ~;90° 
7 0.07 e . 0.02 e@ 
| t.37 26-4 
$4 4.26 +) $+ 4.26 -| 7.37 (1-138) 
OMIA ele > cos ST ation rach) (1-139) 
and the total response is 
~4 26% 
x(t) =-476 e cos (4.474 +2.4 rad) + 
+ 0,/4 ete cos ( 7.47 + 4.57 rad ) (1-146) 


(2) Unit Step Response 
The spring mass damper is driven by the singularity 
input 
FQ) = ult) (1-141) 
The differential equation becomes: 
Mer wre, x =) (1-142) 
Neglecting initial conditions and Laplace transforming, we 


Obtain: 


Xart (s) # (s*+ 25uns tun) = s (1-143) 
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Or 


Xoart (s) = S/s?% 25uns + Wn" 
For our example with 5=1 
8.5 
SS 


(54854)? 6+ 85) C5451)? 


-8.51t -8.5it 
Xpary 1) = € - BSite 


and the total response 1s 


-851¢ -8.511 ~g5it 
u(t) = (4+ 66.044) e +e -t%e 
With 5=1.5 
AF ~ 
Econ Cel = ae 
($+3.26)C5+22.28) $+ 4.26 $+ 22.28 
- 3.26 PL GRE RS 
at) icine > aes = EWE 
and the total response is 
-3.26t ~22.28t ~3.26%t z 
x(t) = 6.36 e -2.36¢@ olf C -1i#e se 
With 5 =9.5 
-30.11° We 
5 0.58 ¢e 0580" 


Xpart ($) = —————— 2 


mms ee 
($44.26 +) ¥.37)(S+4.26-j7.37) «#42645 F3F  S44.26-5 4.34 


wget 
Xpary t) = 11Ge cos (7.37 4 + 0.53 rad) 


and the total response is 
-426t 
x(t) = -%%6e cos (F.3F4 +24 rad) 4 


4.26% 
$116 e cos (7.37 +t +0.53 rad ) 
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(1-144) 


(1-145) 
(1-146) 


(1-147) 


(1-148) 


(1-149) 


(1-159) 


Gat) 
Ei P52) 


(1-153) 





1.6 AXES SYSTEMS 

In describing the motion of our flight vehicle, we will 
make the following two assumptions: 

First, we will establish an inertial frame of reference 
on the earth considering the earth to be fixed in space. [In 
this inertial frame, Newton's laws are valid. 

Second, the vehicle will be assumed to behave like a 
rigid body so that translational motion can be described by 
considering it as a particle with all its mass located at the 
center of mass, and rotational motion can be described by 
considering moments about that center of mass and constant 
moments of inertia. 

To describe the translational and the rotational motion, 
one has to establish a second non—-inertial axes system, fixed 
on the vehicle with origin at its center of mass. In this 
reference frame Newton's laws do not hold, in general, and we 
have to add additional terms to our basic fundamental 
equations. 

It is convenient to use cartesian coordinate systems for 
both reference frames used. 

Let's consider the translational motion first. We will 
assume that Oxyz is the inertial reference frame and O'e’y’2/ 
the non-inertial. As in Figure 1-7 shown, fr and fr are 


position vectors pointing to an object and originating from 
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the origins of the coordinate axes correspondingly. fr, 
denotes the position vector of the non-inertial origin 


measured from the inertial origin. 





Figure 1-7 


Inertial and Non-Inertial Reference Frames 
Translational Motion 


The relation between those three vectors is clearly 
fFf-e- £ +7 
with respect to time, we obtain 


gee Les Cheong) 


By differentiating twice, 


Now Since we have assumed that Newton's laws hold on the 


frame Oxyz , we may write 


a (1-155) 
Or substituting Equation 1-154 
ee) See Oo aa (1-156) 


where mis the mass of the object and F the force applied 


OneoO it. 
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Equation 1-156 merely says that we need the extra term 
ant, to describe the translational motion of the object in the 
non—inertial frame. 

Meee @erhen thewrane Oxy?’ is inertial and Newton's laws 
hold. In this case, Equation 1-156 takes our familiar form 

ai eee Gers) 
which 1S similar to Equation 1-155. 

Let's now consider the rotational case. We will assume 

teaktvel! t % =8 so that 0 is the same as 9, as shown in 


Figure 1-8. 


| Kw > 
> 





Figure 1-8 


Inertial and Non-Inertial Reference Frames 
Rotational Motion 


The position vector r can be expressed in the inertial 


frame as 


=xt+tyfezk (1-158) 


“4 


Or in the non-inertial frame 


re x Leyte 2’k eae ey 
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where (x,y,z) and ( x’, y’,2°) are the corresponding axes 


components of Tr. 


Similarly the vector A can be expressed as 


> 


AzAyt+Ayf +A, k (1-168) 


Or 


=» 


A= Ay i's Ay fe Ag ke (1-161) 
Differentiating Equation 1-169, 1-161 with respect to 


time, we obtain 








dA dAx dAy 4 dAz 4 
ye L + == {| + k 


im ~ ai a aR (1-162) 


and 








My dae ne dAyar . dAe ay 
(= a ame ap (1-16 3) 


where (dA/ gt )” denotes that we are taking the time derivative 


of A on the primed non-inertial system where the unit vectors 


AT amt d 


L ee are constants. 


We are interested now in the time derivative of the 


—_> 


vector A relative to the inertial system, but we would like 
to differentiate Equation 1-161 instead of Equation 1-160, 


obtaining thus: 


dA dAx Al dAy al Az 4 ‘ di’ ‘ dk’ 
me SR Sw An See dy are (1-164) 


at 





A 


Or 


> 
21, 


day! yr dt’ pedi’ ye dk’ 
-(F) LA = te * jane (1-165) 


It can be shown that a more convenient expression for the 


three last terms is given by 


4l 





/ dit’ : aj ” eh ~ 7G 
Ay apo fy ne ote pe = wxA (1-166) 


where w is the angular velocity vector of the primed system 
which is assumed to rotate about an axis having the direction 
of ®, through the origin, and with angular velocity ©, 
relative to the unprimed inertial frame. 
For proof, assume for a moment that A is a vector "at 
rest" on the primed system so that 
(dAzn)’ =0 Giaieg) 
We shall show that 
Jax = xk (1-168) 


Using the time derivative definition, we can write: 


dA , kG+at)-AW t 
cise, Ab (1-169) 


Figure 1-9 shows two succesSive positions of the vector A 
which rotates with angular velocity with respect to the 
inertial frame, so that the tip of the vector traces a circle 
perpendicular to the direction of w. 

To show Equation 1-169, we must show that the directions 
and magnitudes of the vectors on the left and right-hand side 
are the same. 

Notice that the direction of the vector lA /ett , by 
Equation 1-169, is the same as the direction of the vector 
A(4+ At ) = AQ) , Which is the same as the direction of the 


vector wxA by looking at Figure 1-9. 
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y= Kear) - AG) 
A(++4t) 


Figure 1-9 


Time Derivative of a Rotating Vector 


The magnitude of the left-hand side of Equation 1-169 


iss: 
[A (t+ At) - Al) (Asing )(w At) 
ee hin ee 

Or 


A 
2). Awsing on 
which is the same as the magnitude of the right-hand side 
vector, |wx Al. 
Hence, both directions and magnitudes are the same, so 
Equation 1-168 holds and if A is not stationary with respect 


to the non-inertial system, we have according to 


Equation 1-165 


43 





_ sy 7 
cA (| ae 


at Vk (1-172) 
which says that the time derivative of a vector A relative to 
an inertial frame equals the time derivative of the same 
vector measured at the non-inertial frame plus the term A 
where © is the angular velocity of the non-inertial as seen 
by the inertial frame. 


Combining translational and rotational motions, one can 


eaSily arrive at the equation: 


A gpl i 
= (4 ~*~ wxA + %, (1-173) 


where f is the linear velocity of the non-inertial with 
respect to the inertial frame. 
One step further, differentiating once again 


Equation 1-172, with respect to time, we obtain 


4 


d -4(¢ 4 ~ dA dw sae 
mia: > ay =] AT ee 
RY! (BY, al AA) anh) $B 2 
2 ( 5 -ix(S ale +WKA? + x A (1-174) 
or ~~ - - 
eed 1 ef AA’ = f= fy Udo 
ae (HF + 2ix(F) es ce = ae (1-175) 


ool 


which says that the second time derivative of a vector A 
relative to an inertial frame equals the second time 
derivative of the same vector measured at the non-inertial 


frame, plus some extra terms involving ¥,A, %%/dt , and 
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(dA/dt )° where ® is again the angular velocity of the 
non-inertial as seen by the inertial frame. 

Combining translational and rotational motions one can 
easily arrive at the equation 
JPA (dA) (day, = (ao iB Oe 
cH? \ le +2Wx\F] + B x(x KA) ese xA +h (1-176) 


where to is the linear acceleration of the non-inertial with 





respect to the inertial frame. 
A special significance of Equations 1-173 and 1-176 is 
given by substituting A with the position vector Y of an 


object. They are taking now the form 


dr dr loxY a 
Tae ie Glas 
or 
= (V) +@xrer, (1-178) 


a 


where V is the velocity of the object, and 


2-» ee - ae 
<t - ee. + 25x(8) + dx(x7) - SB y eo fe (1-179) 
or 
A <CG) + 20 x(V) + Gx (MxF) + Ger e | (1-189) 
where G@ is the acceleration of the object. 
One can recognize that 
9Hx( Vv) = corriolis acceleration (1-181) 


and exists if the object has a velocity (V)’ with respect to 


the non-inertial frame. Also, 


wx(@xr) = centripedal acceleration (1-182) 
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which is always directed toward the axis of rotation and is 


perpendicular to that axis, as shown in Figure 1-19. 


—~ 
w 
>ox (xr) 
AXIS of wx 
rotation 


Figure 1-19 


Centripetal Acceleration 


Since we have assumed that the unprimed system is 
inertial, we may write 
ma =F (1-183) 
and substituting Equation 1-189 
an (&)! = F -2m[dx(V)'] -m[ dx (Bx7)] -m(OxF) - ma, (1-184) 
where Wis the mass of the object and F the force applied to 
ie. 
Equation 1-184 merely says that we need four extra terms 
to describe the translational and the rotational motion in a 
non-inertial frame. 


It is clear that the forces: 
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2m] wx(V)'| = corriolis force (1-185) 
m | ox (dx7)] = centripedal force (eo } 
m(@xt) = Euler force Sea) 
md. = D' Alembert force (Poss) 


are not true forces. They only exist due to the fact that 
our object 1S moving on a non-inertial frame. 

According to our first assumption that the earth is the 
inertial frame of reference, we consider that it is rotating 
with fixed angular velocity W so that we can write 
Equation 1-184 as: 

Pai. =) E - 2rd) “ia (conan) nae (1-189) 

We have seen, up to now, the effects of considering a 
non-inertial frame of reference and our next step will be to 
orient our flight vehicle in that frame. The reason for 
selecting a non-inertial frame will be greatly appreciated 
when one will see how easy it is to formulate equations of 
motion of a body attached in this frame and moving with it. 

For flight control applications, it is appropriate to 
select the origin of the axis system on the center of mass of 
the vehicle. 

From then on we will orient the nose of the vehicle 
towards the positive x-axis, which will also be a symmetry 
axis. The positive y-axis will point to the right and the 
positive z-axis will point downwards having the direction of 


gravity, if one considers straight level flight as 
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Figure I-11 shows, so that the xz plane will be a plane of 


symmetry. 


Vector quantities associated with this axis system are 


Shown in Figure 1-ll and are described as follows: 





VERTICAL 
y AXIS 
NG ~ 
ATRK S 
M positive pitching th SI wwe 
q pitching velocity oe X> nH? 
Y aerodynamic force 
v velocity 
<S 
CENTER 
_ OF GRAVITY 
W positive yawing moment 
L positive rolling momen? r yawing velocity 
* p rolling velocity Z aerodynamic force 
X aeradynamic force w enenity 
v velocity 2 


Figure 1-ll 


Axes System 


Velocity of the center of mass with respect to the 
inertial frame measured on this axis system 
V sutevfewk (1-196) 
Angular velocity about the center of mass 


® = pt+qjerk (1-191) 
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oem force acting on the center of mass 
F = Se On +Zk ino) 
Moment about the center of mass 
G =Lt+MZ+Nk (1-193) 
Directions as shown in Figure 1-1l 
Also, we can define: (see Figure 1-12) 
Angle of attack: 


= tan (WW) (1-194) 


<4 





trace on x2 / 


plane 


Figure 1-12 


Angle of Attack and Angle of Sideslip 


Angle of sideslip 
=! 
= Oty] ( v 
p Nv) (1-195) 
The reference axis system has been adopted by B. Etkin, 
D. McRuer and others, and is usually referred to as body axis 
System (See references). The notation used is quite common 


On most references. 
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Another axes system adopted by C. Perkins is shown in 
Figure 1-13. It is usually referred to as wind or stability 


axes system. 


x’ 


(bod = 
y axis 
nD ai ‘ 

x Cwind) d ae. 


2 
( body) 


Z Cwind) 


Figure 1-13 


Wind Axes System 


In this system the x-axis points to the opposite 
direction to that of the relative wind. Notice that in the 
undisturbed flight the body axes coincided with the wind 


axes, but in the disturbed flight that may or may not happen. 


1.7 RIGID BODY DYNAMICS 
Rigid body is a system of particles in which the distance 
between any two of them is not changing regardless of the 


acting forces. For our flight vehicle the rigid body 
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assumption implies that as far as its motion is concerned in 


three dimentions, only negligible errors are introduced. 

We will consider first the planar motion of rigid body 
where we will find the kinetic energy and the angular 
momentum expressions. 

Assume that the rigid body in Figure 1-14 can move in the 
xy-plane and can rotate about an axis perpendicular to that 


plane passing through the center of mass c. 





Figure 1-14 


Rigid Body - Planar Motion 


We denote ¥ , the position vector of the center of 
mass c, rf the position vector of any particle on the rigid 
body and t’ the vector from the center of mass to that 
particle. 

Usually, the mass of the rigid body is not uniformally 


distributed over the entire volume V. It is convenient to 


SIL 





think for a moment that the entire volume of the rigid body 
containing a mass ™, is projected on the xy-plane on an 
area A. Let's denote m the mass per unit projected area of 
the body. It is clearly seen that the summation over the 
Siti re area of all those masses mca gives the total mass of 


the body. This mathematically is expressed as: 


ms [ne (1-196) 
From Figure 1-14 we have: 
Coe (1-197) 


For a differential mass m dA at which the r vector points, 

} ° ; 
the kinetic energy is expressed as = mda r? » Integration over 
the entire area A gives the total kinetic energy of the rigid 


body. 


i DF 
it i“ mda (1-198) 


The quantity 7° according to Equation 1-197 is expressed 


as 


4° 


eke = (har )(rer) = n2er eon’ (1-199) 


me | 


Then Equation 1-198 becomes 


3 an ee: alan 
k= hi [adh + 3 ferydd of Rada (1-269) 
A ° : 


The last term is zero since 


[Fema = [Fadl [gdh -%m-fam =0 (1-261) 
A A A 


34 





; * / 
We now express the velocity Yr as a sum of two 


velocities: 


oe (1-202) 
where dr’/dt is the radial velocity, @ is the angular 
velocity of the body, 6 the tangential velocity of the 
particle , ? and v denote the corresponding unit vectors in 


polar coordinates, 


Noticing that dr‘/dt =9, Equation 1-208 is written 


J we ae 2 1,2 
Ke 5m 70° a dA (1-203) 


We define 


2 
I -P mA (1-204) 
the moment of inertia of the body about an axis perpendicular 
to the plane of motion through the center of mass. Then 
mr,” 4 


7 A 2 
K = 19 (1-265) 


NI— 


2 

This equation merely says that, in addition to the term 
5 mr that describes the translational motion, we need the 
extra term 119? which describes the rotational motion of 
the body around an axis perpendicular to the plane of motion 
through the center of mass. 


The angular momentum of the same differential mass dA 


about the origin is expressed as 
FxydAr (1-296) 


i) 





Integration over the entire A gives the total angular 


momentum of the rigid body. 


h = | CRF) qd (1-207) 


The cross product term can be written 


ryPe(per)y(ner) = hue ae xr + run +h xr’ (1-208) 


Then Equation 1-287 becomes 


: iat [ach + [Rat ad view | Fqd tx | tnd (1-209) 


In this expression the integral of the last term is zero 
as Equation 1-291 shows. By a similar argument the third 
term is also zero. Noticing that r=arOd, Equation 1-289 


becomes 


“> 


, AR yt) Jam + 0G 9) )ndA (1-21 9) 


If we use the center of mass at the origin Ct =§) we are 


left with: 


lupe 6 | r(#x3) adh (1-211) 


The cross product in the integral indicates that the 
angular momentum vector is perpendicular to the plane of 
motion. Since (Fx 9l =f, the magnitude of the angular 


momentum is given by: 


h = 6 [end eae (222) 


A 
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So the magnitude of the angular momentum of a rigid body 
rotating in a plane, about an axis, through the center of 
maSS, 1S given by Equation 1-212. The angular momentum 
vector points in the direction of the axis of rotation that 
can be easily found by the right hand rule. 

In summary, we have seen that the motion of a rigid body 
in planar motion can be described by a translational motion 
of the center of mass, plus a rotational motion about an axis 
through the center of mass. 

We will now look at the motion of a rigid body in a three 
dimensional space. It generally has six degrees of freedom 
1f not constrained. Hence, we need six coordinates to 
specify the motion. The three of them can describe the 
translational motion of the center of mass of the body and 
the remaining three, the rotational motion. 

We consider two coordinate systems. An unprimed system, 
as shown in Figure 1-15, whose origin is O, for simplicity is 
fixed on the rigid body but its axes do not rotate with 
respect to an inertial frame, and a primed system whose 
Origin and axes are fixed on the rigid body. 

We also assume that the rigid body has a mass 
distribution per unit volume (density) which, in general, is 
a function of the position vector f to a particle on the 
rigid body having a differential volume dv. The total mass 


1s given by: 


2° 








Figure 1-15 


Rigid Body - Space Rotational Motion 
an =| ptf dy ere) 
where p(t) is the density as a function of f and 4 a small 
differential volume where fr points to. 
For the differential volume dy the angular momentum for 
rotation about the origin of the unprimed system 1s given by 
* x] pcF)dy -F | (1-214) 
Integration over the entire volume gives the total angular 
momentum: a — 
h a Oe eal (1-215) 
Bierce t,T are measured on the unprimed system. 
Let ® the angular velocity of the rigid body about the 


Origin of the unprimed system and r the position vector from 


the primed system. 
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It is more convenient to express ¢’, which can be seen as 
the velocity measured on an inertial frame, as 
Y =f +mxr (1-216) 


Notice that 


sod’ a 
r ee ele 


* pte © GeDLG) 


where Y’ and 3 are unit vectors in polar coordinates of the 
primed system. 
mers clear that dr’ /dt =6 according to rigid body 
assumption, and O=5 because there is no relative rotation of 
the differential volume dy that can be seen by the primed 
system. In other words, Yr is a vector at rest in the primed 
system so that we can write 
somal (1-218) 
Substituting in Equation 1-215), we obtain 
h = if Fx (@xF) p(F) dy (1-219) 
Using now the eae product identity 
PeaGExG) = AC) & = (i. Bc (1-228) 
we can write 
rx (wxr) = MWO-(w-F)F Gi=22m) 
so that Equation 1-215 becomes 
h ={ [rtd - (-#)7] pc?) dy Gla) 
In component form we can write 
he=() [rtw, -(a-F) J pet) dy (1-223) 
where t=l1, 2, or 3 denotes the x, y, or z component 


respectively. The dot product is expressed as: 


a7 





> > oan 


(1-22 4) 
Jel 
That permits us to write 
am 
h, = ( fair -fir, |} w: olf)dy 
a il j1 4; ¢ (1-225) 
where dj =1 if i=j and dij =0 if i#} because 
3 ] 
es: SO dmnoew, 2: 
hn ese (1-226) 
We can now define 
ee an eee er = 
Ii, -\'[r di Alle | (OUD) ral (1-227) 


as the ij th component of the inertia tensor I about the 
Origin of the unprimed system. 


So finally, we express the ith component of the angular 


momentum ass: 


we. fe . 
J=t 
In the expanded matrix form this gives: 
hy Tye xy Tez | | ox 
hy = Ixy Tyy Jy2 - | Wy (1-229) 
he Tx. Jyz_s Lez Wz 


This result is compared with the result of the planar 
motion case that is expressed by Equation 1~212, where the 
inertia scalar term is expressed now as a matrix. 

The components of the inertia tensor are given by 


Equation 1-227 in cartesian coordinates as: 
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|e WC ceeyeey - x? | o(r) dx dydz = 


“ff Gies 2?) OCF ) dxdydz 


Similarly, 


2 


Jyy «fff case o(r) dxdydz 


Tan fff C xeu!) OCF) die dyde 


Ixy = Tyx =-fffxy e(*) dx dydez 
Txz < Jzx -- fff o(*) dxdydz 


Also, 


(1-230) 


(1-231) 


(1-232) 


Cb 74878)) 


(1-234) 


Example: Calculate the inertia tensor of a cube of 


Side a, mass m and uniform density o about one of its corners 


(Figure 1-16). 





Figure 1-16 


Inertia Tensor-Calculation of a Cube 


=)8) 





= off (4224) dudyde = = { al ( y2427) dzdydx = 


3 3 
Sf (yvas © )dyds = a (a+ Fa)dx = 
4 4 5 
m7 a a m 2a 2 2 
Similarly, 
2 2 
Tyy =e. = 3 me (1-236) 
aqa,;-a,a 
am 
lay = -e fff xy dxayee ° i] ] pce ee 
av a ot om S ae 
= ={ | xy a dy dx --S/ XZadxe 
__m a @ = - = ma? 
7 ries aaa (1-237) 
Similarly, 
4 
re = oe m a? (1-238) 


T= ma? V4 Jz HA (1-239) 
-'h -' 2/2, 





The eigenvalues and the eigenvectors of an inertia tensor 
have a very special significance. The eigenvectors determine 
the directions of the so-called principal axes of a rigid 
body that pass through the point where the inertia tensor is 
Calculated (usually the center of mass). The corresponding 
eigenvalues are called the principal moments of inertia. 

It can be shown that if a rigid body rotates about one of 
its principal axes, its angular momentum about the point 
which the inertia tensor is calculated, is parallel to its 
angular velocity. 

For proof, let g € principal axis, and B the 
corresponding principal moment of inertia. The eigenvalue 
problem equation which has to be satisfied is: 

Ip =Bé (1-249) 
Or 
liy Pj = BR (1-241) 

Having rotation about a principal axis simply means that 
w has the direction of the principal axis which is expressed 
as: 

b=cB (1-242) 
where C is a constant. In fact, if B is unit vector, C is 
the magnitude of the angular velocity ». 

Now the ith component of the angular momentum becomes 

h. =i ced i, (cf) = ¢ jij Bi =¢ Bf; = to (Celene) (1-243) 
and finally, 

h; = Bu (1-244) 
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In vector form 


Ao ai (1-245) 
which proves that h is parallel tow. 
It can also be shown that 

iL G Any plane of symmetry of a rigid body is 
perpendicular to a principal axis. 

De Any axis of symmetry of a rigid body is a 
principal axis. 

The density @ may not be constant. For example, if the 
xy plane is a plane of symmetry, that simply means that 
o (-x, 4,2) = 6) Cx, tt a) (1-246) 
and if the z-axis is an axis of symmetry, that simply means 
that 
Ge ae eee me pax, une) (1-247) 

3. Since the inertia tensor of a rigid body is 
symmetric, (i.e. Txy = Iyx etc.), the principal axes are always 
Orthogonal to each other (Important result from Linear 
Algebra). 

4. If two or three eigenvalues of a symmetric tensor 
are the same, there exist two or three orthogonal 
eigenvectors associated with that eigenvalue. 

Example: The inertia tensor of a rigid body has been 
calculated with the reference point at the center of mass and 


has been found to be 


62 





40 20 30 

Find the principal axes and the corresponding principal 
moments of inertia. 

Solution: 

Let 4 denote the eigenvector (principal axis) and A the 
corresponding eigenvalue (principal moment of inertia). 

The following eigenvalue equation has to be satisfied 

Tz = Ae (1-248) 


J] 
from which 


A-\5 O 40 


(A-1)7 =0 > oes for 0 (1-249) 
10 0 A-30 


or 
(A-415)(A- 25) (A-30) + 10[-10 (A-25)] = 0 (1-258) 
> (A-25)| (A-45) CA-30) - 100] =o (1-251) 
=» (A-25)(A-35)CA-10) =O (1-252) 


Hence, the principal moments of inertia are: 
A, = 10 A, = 25 Na Se (1-253) 
To find the eigenvectors: 
For A; =186 Equation 1-249 becomes 
2) 0 40 
oO -§ 0 | 4 =0 (1-25 4) 
AO 0 - 20 


6 3 





which yields 
34x i Ve: =0 ae ae = aga ee 
= 19 Day = 0 = May = 60 (l=25:5) 
10, ~ 204.2 =0 a Mx = Qar2 
The eigenvector that satisfied Equation 1-254 is 
2 
4 
It iS more convenient to express the eigenvector as a 


unit vector 


a/"5 
iene (1-257) 
4/ (5 
Similarly, we find that 
O 
9, =] f assosiated with A, = 25 (1-258) 
0 
and 
AG 
ae 0 assosiated with Az ~= 35 (1-259) 
715 


It is important to note that the above eigenvectors are 
Orthogonal since they are coming from a symmetric tensor. 
Since now the principal axes of a rigid body are 


Orthogonal to one another, or can be chosen orthogonal, it is 
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very convenient to have them coincide with the primed system 
axes, the origin of which can be moved to the origin O of the 
the unprimed system. The importance of this arrangement will 
be shown a little later. 

The kinetic energy of the differential volume dv is 
{ 


5 (f-¥)p(F)dv. Integration over the entire volume gives the 


total kinetic energy: 


K 5 (FF) CF) dv 


(1-260) 
Substituting 
F=ox? (1-261) 
we obtain 
re ae (wxt)(wx¥) p(T) dy (1-262) 


Using Equation 1-224, 1-226 and the vector identity 
(Ape Cc) = (A-c)(8-D) - (A.D) (8-C) (1-26 3) 


we can write 


; a 3 
(Oxr }(@xr) =F (> ij 0) Siw) - (2 r all « (04) 
ED acs 
x >} > [dy reid ay tm (1-26 4) 
Sales 


so that 
ae, ou 
K => “a > { Lr? dy HV | Wj | W (1-265) 
Y 
and finally, 
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(ho) (1-266) 


Notice that if the body is rotating about a principal 
axis Chie ) its kinetic energy is maximum K=>hw. That 
happens when the rigid body rotates about a principal axis. 

This compares with the result of the planar motion case 
that is expressed by Equation 1-205, if we neglect the 
translational motion. 

We are now in a position to formulate the rigid body 
equations of motion. 

The translational motion equation is given by the 


equation of motion of the center of mass of the rigid body 
mr, = F (1-267) 
Or 


qi; | ss 
wn ( + wy V = F (1-268) 


The rotational motion about the point (0) which is 
conveniently chosen as the center of mass (c), as Figure l1-1/ 
shows, is given by 

dh = 
aT ( (1-269) 


<P 


where G is the net external moment about the center of mass 


(c) and h is given by Equation 1-228. 
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Figure 1-17 


Rigid Body ~ Space Motion 


The vector dh /dt 1s measured in the unprimed system, but 

it can be expressed as: 
AR en ea 
& (2) + Wxh 


where (dh/at ’ is the same quantity measured in the primed 


(1-278) 


system, and W is the angular velocity about the center of 


maSSe 


So Equation 1-269 is written in component form as: 


ott 


where t refers to a primed axis now. Or better, uSing 


(se) : (xh); ao (29 aa) 


Equation 1-228 
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sd ) _ 
Zu 4) + (xh), ail: 
Jef 


where Ii; 1S a component of the inertia tensor calculated 


( lie 2) 


along the primed axes. 

Notice that the primed axes unit vectors are vectors at 
rest, in the primed system, and obviously Iij is not changing 
with time. Hence, 

3 
Via sey a + (wxh), =; 
= 

This result is interpreted as follows: Suppose the rigid 


(1—273) 


body rotates about a principal axis with constant angular 
velocity W. Rotation about a principal axis implies a 
parallel to h , so that xh =g. This simply means that the 
moment required for this rotation is zero. 

This compares with the planar motion case result which 
comes by taking the time derivative of Equation 1-212. 

eta ae (1-274) 

Equation 1-273 is the equation of rotational motion in 
the general case for rotation about a random axis, where 
external moment is required even for constant angular 
velocity. 

It is now clear why we formulate the equations of motion 
in the primed, non-inertial coordinate system. We have 
avoided the time differentiation of the inertia tensor in 


Equation 1-272 by paying a little for the extra term wxh, 
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In summary, we found the equations of motion of the rigid 
body described in the convenient form which will be used to 


derive the flight vehicle's equations of motion in Chapter 2. 


They are: 
Vv i) .# 
m ( + wxV - F ALaay3) 
and 
3 
———e dw: —~ 
x aes ~ ne 
= qi 7 ‘ (axh)-; = G; (1-276) 
Je! 
where 
=e 
oe | re 
hi = L Tij lon) 
Js= 


and Ij; is calculated about the center of mass with respect 
to the primed axes. 

Finally, if the primed axes are principal axes of the 
rigid body it is not difficult to see that the inertia tensor 
has only diagonal terms representing the principal moments of 
inertia (eigenvalues). This result comes from the 
diagonalization process of a square matrix. So, in this case 


Equation 1-276 is simplified as: 


dw; + 7). 
qi ae (Oxh)i =, (1-278) 
where 
hi = lit %i (1-279) 
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CHAPTER 2 : 


VEHICLE'S EQUATIONS OF MOTION 


2.1 INTRODUCTION 

With the methods of obtaining rigid body equations of 
motion, established in the preceding chapter, we now proceed 
in this chapter to formulate the equations of motion in a 
most generally useful form. 

The equations are basically simple in the sense that they 
are derived by making use of the basic Newtonian laws for 
Forces and moments. Great advantages are obtained, as we 
have seen, by using a non-inertial frame of reference, 
because with an axis system fixed in the body, the inertial 
properties are constant and that simplifies the moment 
equations. 

Some penalties are obtained by introducing non-inertial 
frames, the most important of which is the addition of extra 
terms in the basic equations of forces and moments, in trying 
to make the basic laws valid in this environment. 

In any case, although the equations are simple, they are 
highly non-linear and not generally solvable mathematically. 
The formulation is greatly simplified by the use of 


linearization techniques, and small pertubation theory. 
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This method assumes an equilibrium position which is, in 
fact, a solution to these equations. That allows for the 
changes in forces and moments to be expressed as linear 
Functions, in terms of small pertubations in linear, angular 
velocities, and other parameters involved. 

The equations of motion will be derived for the body axis 
system established in Figure 1-1l and will be left in their 
dimentional form. Non-dimensional form will be mentioned as 
it applies to body axis system (Etkin's) and wind axis system 
(Perkin's). 

Finally, the equations of motion will be split into two 
Separate groups, the longitudinal and the lateral groups of 


equations. 


2.2 FORCE EQUATIONS 
Starting from the vector force equation 


~ 


> _m( 2 «xi 


alt (251), 
where 
Titan + ZY? +IZk (22) 
V = vt ah (2-3) 
b= pirqjerk (2-4) 


we can derive the force equations in component form. 
The cross product term is expanded as: 


ixV = (qw-vr) % + (-pweor) 7 + (pv-qu)k (2-5) 
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and the force equations of motion are Simply 


2x = m Cu + qw-vr} (2-6) 
LY = m(v-pwtur) (2-7) 
XZ = m(w:pv-qv) (2-8) 


These equations describe the translational motion of the 
center of mass. 

Consider the motion of the center of mass in the xz plane 
as Figure 2-1 shows, in two successive different locations on 


the flight path during time ot. 





Figure 2-1 


Motion of the Center of Mass in the xz-Plane 


The acceleration along the x-axis 1S expressed as: 


; (y+du) cos 58 + Cweow) sin d@ -v 
— 3t (2-9) 





For small pitch angle 9, cos0?=1 and sini#®=0, so that 


U+OU + wood + OwdG -y 
OS Ts ci we (2-18) 


qZ 





letting now ét-0 and neglecting the higher order term dwoO, 
we obtain 
ay = V+ wa (2-11) 
and for the corresponding force equation 
2X = m Cd+wg) (2-12) 
Which is the same as Equation 2-6 by neglecting the term vr. 
Similarly, the acceleration along the z-axis is expressed 
as: 


(w+dw) cos30 - (u+dv) sin 60 — w 


Gl. = 
. <I (2-13) 


which is simplified as: ~ 


w+dw - vod - avd9 - w 


a aaa aT aa (2-14) 
Or 
a, = W- QD (2-15) 


and the corresponding force equation becomes: 
le m ( w -qv) (2-16) 
which is the same as Equation 2-8 by neglecting the term py. 
If we now consider the motion of the center of mass in 
the xy plane as Figure 2-2 shows, in two successive different 
locations on the flight path during time Ot , we can derive 


the acceleration along the y-axis as: 


Ay = 
J Sk 


(2-17) 


Us: 








Figure 2-2 


Motion of the Center of Mass in the xy-Plane 


For small yaw angle yp, cosoy=1 and sindw=dv, so that 


(vedv + VOW +Ov0y -v) 


ae St (2-18) 


letting now dl>O and neglecting the higher order term Ovew + 


we obtain 
dy = V+or (2-19) 
and for the corresponding force equation 
ZY = m(V+tor) (2-20) 
Which is the same as Equation 2-7 by neglecting the term pw. 


In summary, the force equations are simplified neglecting 


higher order terms as: 


pak = m (d+qw) (2~21) 
LY = m (yeur) (2=22) 
Dales Ere (w - qu) (2=23)) 
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2-3 VEHICLE SYMMETRIES 

Up to this point we have made two assumptions. The first 
assumption was that earth is considered to be fixed in space 
so that we can establish an inertial reference frame on it. 
The second assumption was that the vehicle is considered to 
be a rigid body. 

Facts are simplified further by considering the 
symmetries of the vehicle. For an airplane or a missile, the 
xZ~plane is a plane of symmetry, if one assumes that the 
following holds for the density distribution. 

o(x,y,2) = @ (%-y,2) (2-24) 
so that 


Ixy =Jyy = ff e (x,y,2) xy dx dydz --|ffels-y,2)%9 dxdy cle 


Notice that Equation 2-24 says that ( is an even function 


(22 51) 


with respect to y but the integrant function is odd because 
it is a product of an even and odd function. Since the y 
integration has symmetric limits, the integral is zero. 


For the same reason 


Tyz=ley--ff e(x,y2) 42 dxdy dz = Q (2-26) 

The inertia tensor, therefore, for an airplane becomes: 
Tx 0 -Iyz 

0 Tyy ) (2-27) 


~Ixz Q J22 


iS 





Notice also that the y axis is an eigenvector of the 
inertia tensor (Equation 2-27) because it satisfies the 


eigenvalue matrix equation: 


lhe QO -l3 O 0 
0 Iyw o 1} = Ty] 1 (2-28) 
-Tyz 0 Tez O 0 


Hence, the y~axis is a principal axis, with Jyy being the 
corresponding eigenvalue or principal moment of inertia. 
This was expected since y-axis is perpendicular to the plane 
of symmetry xz. 

Unfortunately, this is the only symmetry that exists 
generally in any conventional aircraft. 

In addition to this symmetry another symmetry can be 
found in gun projectiles and in most missiles and 
Space-crafts. They have a second plane of symmetry, the xy 
plane, so that one can assume 

p(x,y,z) = p (x, y,-2) (2-29) 
and hence, 
Tyz = Izy =O (2-38) 

Then, the inertia tensor for a two-symmetry vehicle 
becomes 

Tyx 0 O 


O 0 J22 


qs 





Notice that in addition to the y-axis being a principal 
axis, the z-axis, as well as the x-axis, are all principal 
axes with corresponding principal moments of inertia Izz, and 


Ixx , respectively. 


2-4 MOMENT EQUATIONS 


Starting from the vector moment equation: 


* dw > | 
é =| a ye (aoxh) QOS 
Or in component form 
3 
- dw; - Co 
724 mo + (oxh)i (2233 
where 
h, = ij 0; (2-34) 
G eLt+Maent (2-35) 
ho = hy t +hy Ze hk (2-36) 
aw ° A e o Ai 
—- = L 4 L 
ee Preeqier (2-37) 
Equation 2-34 is expanded as: 
hy = p lx mill I yz (2-38) 
he = -p lxz +r zz (2-49) 


also, the cross product terms is Equation 2-33 is expanded as 


(xh) =(qh,-rhy )4 + (rh, -ph; ) 4 + (phy -4 hx) (2-41) 


Dyk 





Or 
(axh) = (-9p Ixz + qt lz - rq i Ls 
7 (rp IR - r? os ae p° ies ~ rp Tes) { 


+ (q p Iyy ma? Txx erg) iat 


(2-42) 
Or 
(@xh) = | -qoTye + ar (Ter-Tyy)] t 
t ip (Txx-I,2 ) t lye (p2-r2) | i 
+ [ pq (Syy-T) 4 rq Ine] & 
L pq, (Syy-To) + rg Tee | = 
Thus, the moment equations are: 
= xx p = Tx2 " - 4 Tye + qr (T22-Tyy) (2-44) 
M = Tyy 4 $ rp ( Ixx-Izz) + Ixz (p? -r?) (2=45 ) 
N = Ia ¢ - Ixz p + Pg (Tyy -Tyx } + FQ Ix2 (2-46 ) 


2.5 EXPANSION OF GRAVITY FORCES, THE LINEAR AND ANGULAR 
VELOCITY 
The total force acting on a flight vehicle can be split 
into a summation of aerodynamic forces (lift, drag and 


thrust) and a body force which is the gravity force. 
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Figure 2-3 shows an accidental orientation of the vehicle 
> 
with respect to the vertical. The gravity force W can be 


expanded in components along the body axes as follows: 


Wy = - mq sin 9 (2-47) 
Wy = mga cos 9 sino (2-48) 
W, = mg cos8 cosh (2-49) 






passes) 
} o 





Y 
8 
ie: ang cos @ 
z 
vertical plane containing x- axis view om yz plane 
and weight vector 
Figure 2-3 


Expansion of the Gravity Force 


For the expansion of the linear velocity vector along the 
body axes consider in Figure 2~4 the inertial reference frame 
being Ox.y,2, . Transformation to the body axes coordinate 
System is obtrained by rotating the reference frame in an 


Order-dependent manner as follows: 


(a) Rotate around reference axis Oz, by the yaw angle », 
which defines the intermediate reference system Oia = 


Ox, yy 2, ° 
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(b) Rotate around intermediate axis Oy, by the pitch 
angle 0, which defines a second intermediate system Oxsy525 = 
Oxy,Z5 « 

(c) Rotate around final axis Ox by the roll angle $, to 


define the desired body axes system Oxyz = Ox,y2 , 





8 > pitch 
b -» roll 
yr yaw 





Figure 2-4 


Expansion of the Linear Velocity 


The components of the velocity vector in the body axes 
System (U,V,W) are defined in terms of the vector components 
in the fixed inertial frame (U,,V, ,W,) by means of the angles 
9,6,¥Y and the intermediate vector components (U.,V,,W,) and 


(U;,V3,W3). 
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The following relations hold: 


for the first rotation 


UV, = U, cosy + Vi siny (2-50) 

Vv, = -U, sinw + V, cosp (2-51) 

W, = Wy (2-52) 
for the second rotation 

U;= UV, cos 8 - W, sin (2-53) 

V3= V2 | (2-54) 

W3> U, sm@ + W, cos 4 (O-8S > 


and for the third rotation 
UV =U, (2-56) 
Y = V3 cosh + Ws sing (2-57) 
W=-Ys sind +W3 cosd (2-58) 
so that the velocity components along the body axes are 
expressed as 
U= UV, cospcosO + Vy sinweoos® - Wi sin (2-5) 
Y=, (cosy sin® sing - sny cs 6) +V, (snp sinO sind + cos pcos d ) 
+W, cos@ snp (2-68) 
W =U, C cosy sin® cosh esmysmd) +V, (simp sin cosb -cosp sind ) 
+ W, cosO cosh (2-61) 
Similar approach can lead us from the body axes system to 
the inertial for which 


UV, = Voosd cCOSw + VisindsmO cosp - cosh siny ) 
+ W( cos sing cosw + sind sinp ) (2-62) 
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Vv, = VeosOsiny + VC sng@smOsiny + cosh cosy) 
+ W (cosh sin8 smMp - snd cos w) (2-63) 
W, = -Usin®8 + V smdos@ + W wsdwsd (2-64) 

Integration of Equation 2-62, 2-63 and 2-64 give the 
expressions for the translational motion but in other than 
Special simple cases they are not integrable. 

Following a Similar procedure we can expand the angular 
velocities p,;q,r into components along the body axes. 

With the vehicle in horizontal flight, one can see the 
yawing-velocity vector having the same direction as gravity, 
the rolling velocity vector having the direction of the 
longitudinal axis and the pitching velocity vector, 
perpendicular to both as shown in Figure 2-5. Those 


Orientations are obtained according to right-hand rule. 





2 
vertical plane reer view om y2 plane 
and weight y 


Figure 2-5 


Expansion of the Angular Velocity 
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Then the angular velocity vector wW can be expanded in the 


body axes as follows: 


P=d-wsind (2-65) 
oF 8 cosh + wcosO sIn® (2-66) 
= -§ sinh + p cosO cos (2-67) 


One can now solve the Equations 2-65, 2-66, 2-67 for om 
0,9. 
Multiplying Equation 2-66 by cos$¢, Equation 2-67 by 
Sind, and subtracting, we obtain: 
Q = q 60s > -r sing (2-68) 
Multiplying Equation 2-66 by sing, Equation 2-67 by 


cosg, and adding, we obtain 


» sin® . cos 
- 3 Geel | | gee (2-69) 


Then substituting Equation 2-69 into Equation 2-65, we 








obtain 
b =Ppt+q sind tan® + rcosptanO (2-79) 


Integration of Equations 2-68, 2-69 and 2-79 give the 
expressions for the rotational motion but in other than 


Special simple cases they are not integrable. 


We can now modify the vehicle's equations of motion by 


considering the gravity as follows: 


Force equations: 


Xx =m ( v +qw -vr +gsin® ) (2-71) 
Y= m(v +ur-pw - gos 8 sing) (2-72) 
Z= m(w+opv “qv -3 «58 ws) (2-73) 
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where now the X,Y,Z forces represent only aerodynamic forces 
(i.e., contributions of lift, drag and thrust). 


Moment equations: 


L=Iyxp -Ixer - gp lx + qr (Tez -Iyy) (2-74) 
M = Iyy q +9p (Ixx ~To | + Ly Cpr =r) (2-75) 
N= Iz 7 - Ixz Pp + py CTyy-1xx ) x rq, ixz (2-76) 
with auxiliary relations: 
Dp) b - psin (2-77) 
q = 8 cosh + wcos9 sind (2-78) 
r=-Osnd + w csO sd eso) 


Notice that the gravity force contributes zero moments 


because it acts at the center of mass. 


2.6 LINEARIZATION OF FLIGHT VEHICLE'S EQUATIONS 


The three force and the threee moment equations derived 
in Section 2.5 must be linearized for use in the stability 
and control analysis. It is assumed that the motion of the 
vehicle consists of small deviations about a reference steady 
fmenant path. 

A steady flight path is described by the general 
equations of motion, if we consider all time derivatives 
zero. We denote by capital letters and the subscript (0), 


the steady flight conditions, obtaining thus: 


X = m(QoWe - VoR + g sin Oo) (2-88) 
Yo = m(U,R, - BW - g cos © sin ® ) (2-81) 
Z, = m (YoYo — Ges - g cosOo cos &, ) (2-82) 
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Lo = ~Gof Ixz + @QRo CIz2-Iyy) (2-83) 
Mo Cae ( Txx-I,, ) + Iyz GPs =k. ) (2-84) 
No = VoQo CIyy-Txx) + 8.Q Ixz (2-85 ) 


By subtracting from the original equations the above 
Steady flight equations, we obtain the perturbed equations. 


For example, subtracting Equation 2-88 from Equation 2-/1 we 


gets 
X-Xo = am | U+qw - GoW -vr +VoRe + g ( sin@ - sin Oo) | (2-86) 
Notice that 
qw = (Gotdq )(Wo tow) = QoWs + Wo dg + Qe dw rs bese 
(2-87) 
Similarly 
vr = VoRe + Vo dr + Rody (2-88) 
Also 
sin = sin (@,+68) = sin@ cs 68 +cos®@ sin (2-89) 


For small perturbation quantities, Equation 2-89 becomes: 
snQ = sin@ + (cos) dO (2-90) 
Combining all these, Equation 2-86 becomes 
oS am | D+ Wo dq, +Qoow -Vo dr - Ro dv + (qcos Oo) 58 | (2-91) 
Although we are deriving perturbed equations, it is more 
convenient to drop the differential symbol Q. 
Similar work leads us then to the following set of 


linearized equations: 
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he 


* 


X= m[v +Weq +Qow -Vor - Rov + ( g cos Oo) ©] (2-92) 
Yeum[v+ Urt+ Rov -Wop-Rw - (9 cos Q cos &,) & +(gsin@ sini, 9} (2-93) 
Z=m | w + Vop ria a8 - Veg -Qov +Cq cos OQ cos & ) d +(g sin ® cos dell 2-94) 


L=Ixxp -Ixz fF - (Qop +feq)Ixz t Gor + Rog )(Tzz-Tyy ) (2-95) 
Ve Tyy q+ (Ro p+ ter) (Ivx-Iez ) +Ixz( 2Pp - 2Ror) (2-96) 
N = Jazr - Ixz p + {PR a+Q.p)CIyy “Txx) + (Rog +@or) Iyz (2-97) 


In addition, for the auxiliary relations described by 
Equations 2-77, 2-78, 2-79, the corresponding perturbed 


equations are: 


p= 6-wpsin -(Y, cos O,)O (2-98) 
q = Q cos b, -(R sin@, sm ®)O +(W cos ©, cosh, - O, sin) b + 
+ w cos ©, sin Do (2-99) 


r= pcos 0, cos ® - CW, cos sin dy +, cos Do)  - Osin & 
- (0, sm®, cos ®) 8 (2-166) 
Although we have obtained linear equations, they are 
extremely complex and not suitable for control analysis. 
A further simplification is obtained if we consider the 
reference flight path to be straight, l.e. 
V, = 0, =O (2-101) 
with no angular velocities 
Wy = Gly 2 Ny =e (2-162) 
Wings level flight, i.e. 


, = O (2-103) 
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symmetric, i.e. 
Wo =O (2-194) 
Vo =O (2-185) 
With this assumption, the reference flight path is 


described by the equations: 


%, = mag sin Sr (2-196) 
\oueao (2-187) 
Zo = - mq cos O (2-188) 
Lo = Mo = No =0 (2-199) 


while the perturbed equations become 


X=m]d+WMq + q cos ©, 2 | (2-119) 
Y=m|v+Uor+Wop -g ws 4] (2-111) 
Z=m|w-og +t gsinOo9 | (2-112) 
L=Ivxp-Ilyr (2-113) 
M= Tyy 9 (2-114) 
N=, ¢ - le 4 (2-115) 
with auxilliary relations: 
p= >-psin Os (2=116) 
the 9 (2-117) 
r= Y cos 0, (2-118) 


We find that the equations are separable into two groups. 
The first group associates the equations that describe the 
motion in the xz body plane, called longitudinal equations 
and the second group associates the equations that describe 
the motion in the xy body plane, called lateral directional 


equations. 87 





We find also, the force equations divided by mass and the 
moment equations divided by the corresponding moment of 
inertia. 

In the longitudinal equations the independent variables 


ares 


Vv x-axis linear velocity 
W = z-axis linear velocity 
q = pitching angular velocity 
In the lateral-directional equations the independent 
variables are: 
V = y-axis linear velocity 
p = rolling angular velocity 
r = yawing angular velocity 
With all these quantities being perturbations from the 
initial flight conditions. 


We arrive then to the following form of equations: 


Longitudinal equations 


| : 

ae = D+ Weg + gcos & © (2-119) 
4 . We 

ee ee (2-128) 
| : 

cm 4 aie 


dyy 


Lateral-directional equations 


Y =v +Ur-Wop -gcos & (2-122) 


i 
m 
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a, - 2 ; 

ies as Co (2-123) 
{ tyz 

S22 N es P (2-124) 


For further simplification, it is convenient to introduce 
here the wind axes system, which as we have seen is produced 
by a rotation about the y-body axis so that the x-wind axis 
points in the opposite direction of the relative wind. This 
axlsS 1S always tangent to the flight path and makes an angle 
meme tignt path y, with the horizontal, as shown in 


Figure 2-6. 





Figure 2-6 


Direction of the Body Axes with Respect to Wind Axes in 
Straight Level Perturbed Condition 


Since 


Q=yto (2-125) 
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in the steady situation 
Oo = Yo + We (2-126) 
where 
Oe (2-127) 

In fact, in order to maintain straight level flight, % 
has to be small. The appropriate relation for expressing the 
angle of attack « in terms of the w-velocity if obtained by 
Equation 1-194 if one considers small « and w, so that one 
can express: 

a = W/U, (2-128) 

From Equation 2-127 then we deduce: 

a, 2 Wey, 20 > We 2 (2-129) 

Notice that by this alignment, since the wind axes are 
Fixed for fixed relative wind, we can measure the 
perturbation quantities about y, in the wind frame. 

It is further desirable to express the Lateral 
Directional equations in terms of the sideship angle 6 
instead of the v- velocity. The appropriate relation for 
saoeessing the vangle 6 in terms of the v- velocity is 
obtained by Equation 1-195, if one considers small 6 and v, 
so that one can express 

B= Y/U, (2-138) 

So introducing 

from Equations 2-126 and 2-127 
Qo = fo (2-131) 


90 





from Equation 2-129 

W, =0 (2-132) 
from Equation 2-136 

v= BU, (2-133) 
we obtain 


Longitudinal Equations 


— = D + g cosy, © 

—Z = w-Uq+gsiny.O (2-134) 
=i" ene 

quy $ 


Lateral Directional Equations 


Y=v+uUr = J COS fe b 


4 
am 








ep a 2-135 
Txx - Ixx ( 
1 Ixz 
= ( --- 
J22 _ slaps 
with auxilliary relations: 
Ds ob - y si” Ye 
q = Q (2-136) 
; oe p COs y. 


The only thing that remains is to express the left~hand 
Side forces and moments in series of terms involving 
aerodynamic derivatives and small perturbation quantities, 


that in fact generate the forces and moments. 
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2./ EXPANSION OF FORCES AND MOMENTS 

In Section 2.6, it was assumed that the motion of the 
vehicle consists of small deviations around a reference 
Steady flight path. 

Small disturbance theory applied on the perturbed 
equations of motion assumes that the aerodynamic forces and 
Moments are functions of the disturbance velocities (angular 
and linear), the control surface deflections and their 
derivatives. 

In general, for a function 

OARS... ) (2-137) 
the Taylor series expansion about a reference value 
Sol dbs 5 (2-138) 
assuming that the function has derivatives of all orders 


defined, is written as 
K = re) = cs ok 
Ko +3K ke + (55), 84 + (55 


fe 
“(CB Om oa Pe) 95 (FE) owe]. 


where the derivatives must be evaluated at the reference 





(2-139) 


value denoted by the subscript (0). 
Now, making the assumption that the deviations 64, dp,a 
dv 


, are so small that product terms are essentially zero, we 


obtain 


a 


SK = ) 3 (E) 5p +( PE) Sve 


OA Ov |, (2-148) 
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Using this method and provided that the above assumption 
is valid, we succeed to express the perturbation quantity 4é« 
as a linear combination of the perturbation parameters from 
which it depends upon. 

People in the field like to denote the partial 
derivatives as Ky, Ky, Ky, etc., instead of (0K/a1)., 
(OF J/op ), (OK /Qv ), , etc. In addition, according to the 
convention met in Section 2.6 of dropping the 5, we obtain: 

K=K,A+ Kp p+ Ky vt. (2-141) 
where Kar Ku, K,, are called stability derivatives. 

Fortunately, not all the stability derivatives in an 
expansion of a force or a moment have a value. Some of them 
are zero or essentially zero. 

Expansion then of forces and moments, in terms of the 
Significant stability derivatives yields: 


For longitudinal equations: 


=X = Xu + Xyw +X58 (2-142) 
— 7] = ena +/,W +Zy W + Lo 4 +736 (2-143) 
i ; | 

jig tt ea la + Ma 4 + Mso (2-144) 


where it is clear that 


Cian (2-145) 





and denotes the control input due to control surfaces 
deflections. 


For the lateral-directional equations: 


—Y =v +i +YpP +Yrr +Y58 (2-146) 
FL =leb+lae rlpp tlrr «lsd Nan) 
= N = Ne@ + Nag +Npp +Nrv + Nod (eae) 
where it is clear that 
Yy = a , Yee = (2-149) 


Therefore the equations of motion become: 
Longitudinal Equations 
V+ Xpv-Xyw + g cos ¥, © = X50 
-Zyu +0 -2Zyw-Z,w -Uoq -Zqq +gsnge = 258 (2-158) 
- My v - My w —— My w 4g - Ma 4 = Ms 9 


Lateral Directional Equations 


v-Wv-Yyi -Ypp + Ur -Yer -gcosy.g = Ygd 
-lep - Le 6 +p -Lop -Ia, ¢ -ler = 158 
-Neé - Nee - Iigq - Nop +r-Nrr = Nod (esti ily 
With auxilliary relations of Equation 2-136. 
It is more convenient to Laplace Transform the equations 


and arrange them in matrix form. Further augmenting the 
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equations by the auxilliary relations, using Equation 2-138 
for the lateral directional set and dividing the Y-equation 
by U we obtain the form 


Longitudinal Equations 





5-Xy =u 0 J COS Xo v X5 
-Ly s-slw -Lw -(Up+ SIN ¥o Ww bs; 
“i gsm : S (Ga132) 
~My -sMw -Mw s-My O 4 M§ 
l 0 0 l => | | 8 0 
Lateral-Directional Equations 
Pr « g + 
$- s Yy -Yv -Yp 4-Y, = v. COS Xo O B Y5 
aj ] 
- slg -lg s-Lp -$ = at O O P L§ 
“sNi-Ne 5 BE-Np Ny 0 re Pane poms 
O 4 0 = s(sing.) | | 0 
0 0 t O -$ (cos Yo) : q 9) | 
where 
Yo | Yp LU, (2-15 4) 
=a Gai (2-155) 
an (ee (2-156) 


We will further simplify the equations of motion in the 
following chapters in a form that is more suitable for 
obtaining aerodynamic transfer functions, and computer 
solutions. 

Our interest first will be focused on examining the 
aerodynamic stability derivatives. They will be left in 


their dimensional form rather than non-dimensional. 
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CHAPTER 3 


3.1 INTRODUCTION 

Stability is a very important concept in successful 
vehicle flight. We have encountered the terms "stability" 
and "equilibrium" but we did not emphasize particularly those 
concepts. 

In the last chapter we mathematically obtained equations 
of motion and linearized them about a referenced straight 
flight which was assumed to be an equilibrium condition. 

But what is equilibrium? Equilibrium is a state of rest 
Or uniform (unaccelerated motion) in a straight line where 
all forces and moments acting on a body are balanced out. 

Stability 1s a property of an equilibrium state. We will 
consider two types of stability which are of great interest 
in flight vehicle dynamics, defined as follows: Static 
Stability is defined by the direction and magnitude of the 
initial tendency to return to equilibrium condition after a 
small disturbance has occurred. It is distinguished by 
positive static stability if the body tends to return to its 
equilibrium state, negative static stability if the body 
diverges from its equilibrium state, and neutral static 


Stability if the body remains in the disturbed position. 
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As an example, the pendulum considered in Section 1.3 has 
a poSitive static stability at its Q=0° position but has a 
negative static stability at its 9=188° position, while its 
motion remains invariant if the pivot is moved in any other 
position and is said to have a neutral static stability. 

Dynamic stability determines the resulting motion in time 
of the body if initially disturbed. It is distinguished 
Similarly in positive, negative, and neutral dynamic 
Stability. 

Various conditions of dynamic behavior were considered in 
Chapter 1, in analyzing spring mass damper systems. 
Figure 3-1 summarizes all possible time responses of any 
dynamical system. 

Dynamic stability implies the existence of static 
Stability, but the opposite is not necessarily true. 

Now the geometry of the vehicle comes to play a unique 
and important role in specifying any particular mode of 
dynamical motion. All surfaces exposed to the airstream 
contribute forces and moments about the center of gravity, 
which are expressed as we have seen in the last chapter as 
linear combinations of stability derivatives. 

The physical meaning of each of the stability derivatives 
is straight forward. Thus, X, is the rate of chance of the 
force generated in the x-direction due to a v-velocity 


perturbation, divided by the mass of the vehicle, while Xu v 


oF) 
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Figure 3-l 
Various Responses of a Dynamical System to a 
Small Disturbance 
a,b,c Non-Oscillatory Modes, d,e,-f Oscillatory Modes 
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is the force generated in the x-direction due to a y-velocity 
perturbation divided by the mass of the vehicle. 

In this chapter we will give more insight in evaluating 
the stability derivatives in terms of the vehicle's geometry 


and airstream properties. 


3.2 LONGITUDINAL STABILITY DERIVATIVES 

Longitudinal stability derivatives are encountered in the 
longitudinal equations and, hance, they uniquely specify the 
X-force, the Z-force and the M-pitching moment. 

Dimensional analysis approach indicates that we can 
express the aerodynamic forces and moments in terms of their 
corresponding non-dimensional coefficients which are usually 


deduced from wind tunnel tests of models. 





Figure 3-2 


Forces and Moments Considered in Longitudinal Motion 


a2 





In particular, we can express the longitudinal forces and 


moments as 


ites L=C. 5 VS 62 

drag D=G5eV'S (Geo) 
! 2 

pitching moment M=Cm poe V"s8¢ (3-3) 


where C,,C,,Cy are the corresponding lift, drag, and pitching 
moment coefficients, p the air density, V the velocity of the 
vehicle, S is a reference area uSually taken the wing area, 
and c, a reference length, usually taken the length of the 
mean aerodynamic chord. 

The Longitudinal Stability Derivatives can be evaluated 
as follows: 

1. The Xv Derivative 

Is the rate of change in the X-force due to a 


UV=-velocity perturbation divided by the mass m. 


t+ OX 
Ay = = (3-4) 


am Sv. 
According to Figure 3-2, the X-force is expressed as: 
X=Lsna -D cosa + ' cos x, o=)) 
where T denotes the thrust. 


For « small, “«, essentially zero, and making use of 


Equation 3-2 we obtain: 
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| Z 
se eS +1 (3-6) 


Notice that because a small change only in v-velocity 
occurs 
V> = (Uo+u) + v4 w? (3-7) 
and since products of perturbations are eSsentially zero 
V2 = Uy + 2Uou (3-8) 
where Use is the reference flight forward velocity which for 
practical purposes is assumed equal to V. 
Then Equation 3-6 becomes 
X=-5 0 (Ue +2) SC, +1 ee 
Neglecting the effect of the v-velocity on the thrust 
T, by differentiating Equation 3-9 with respect to uv, and 


dividing by the mass m, we obtain 








mo > om Sy Bid) 


mM 2 v (3-11) 
Finally, if we define 
Dee 
, ~ 2 vu (3-12) 


we arrive at the following final form of the Xy derivative 


OSU, 
leer (Co +O,) (ie) 


Orr 





The quantity Cp, is essentially zero for low subsonic 
speed, but sometimes near the critical mach number 
(9.8<M<1.0@) where a large increase in drag occurs, Co, 
obtains a large positive value. 

2. The Zy Derivative 

It is the rate of change in the Z-force, due to a 

V-velocity perturbation divided by the mass m. 


LWA 
Ov (3-14) 


= 4 = 


Ly = 


According to Figure 3-2, the Z-force is expressed as: 
Z = Woos@ - Loosw -Dsina - Tsina, (2s) 
which is simplified to: 
Z=Wcos® -L (3 =116)) 
Neither W or @ depend on V-velocity. Since by 


Similar approach as in the X»y derivative, we obtain: 





Ly Boe (C4. +C.,) (3-17) 


where 

(3=ak3)) 
The quantity (C, is essentially zero for low subsonic 

speed, but sometimes near the critical mach number may reach 

a considerable value and its sign can be changed depending on 

the airframe geometry, mach number, dynamic pressure and 

aeroelastic properties. 


For subsonic straight level flight (¥y=0) 


LO2 





aie ee 
"pve $/2 (3-19) 


and Z is expressed simply as 


__ 24g 
fy = - a (32005 


3. The M, Derivative 
It is the rate of change in the pitching moment M due 
to a V-velocity perturbation, divided by the Iyy moment of 
inertia. 


{ OM 
Tyy Ov (3-21) 


The Mv derivative is evaluated simply by 


differentiating Equation 3-3 with respect to v. 


7 9 SUs¢ 
Pris Tyy (Cur + Coa) (3-22) 
where 
G = Up “Cu 
Mo = > (ges 


The quantity °M) is very sensitive upon the variation 
of the v-velocity and its sign can change depending upon the 
geometry, aeroelastic properties, mach number and dynamic 
pressure. Thrust can also affect Cu, for propeller driven 
eubGe@iedit 

4. The Xw Derivative 
Pre rseene rate of change in X—-force, due to a 


w-velocity perturbation divided by the mass nm. 


OS 





Yes 


= oo 


| 
z | $< 


(3-24) 
A w-velocity perturbation causes a small change in 
the angle of attack «. Thus, in this case we cannot neglect 


the angle of attack w which is approximately given by: 


Up (3-25) 
Therefore, the X-force is expressed as: 
X=Ll sine -Deosxe +1 (3=26) 


Or 


oe oV'S(G sma - G cos a) +] (3-27) 


Neglecting the effect of the w-velocity on the thrust 
T, by differentiating Equation 3-27 with respect to w and 


dividing by the mass m, we obtain 





ys 
Os = aus Se = = Si) & -@, cos.) (a8) 
Now aSSuming « small we express 
fe 
(C. - Co, (3-29) 
where 
on _ © 
on (3-30) 


The quantity G,. represents the change in drag 
coefficient due to angle of attack changes. When the angle 


of attack increases, the total drag increases and this 
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quantity has normally a positive value. It is affected 
primarily by the wings, and a little by the tail and the 
fuselage. In aircraft dynamics it is not very important 
parameter mainly because the changes in angle of attack are 
small and, hence, its value is very small. 
5. The Zw Derivative 
It is the rate of change in Z-force due to a 


w-velocity perturbation divided by the mass m. 


Zy == = (3-31) 
The Z-force simillarly is 
Z=Wos8-Leoso -D sino (3-32) 
or 
Z =WosQ - 5 pv's (CG, cosx -G sine) (3-33) 


Neither W or 9 depend on w-velocity, since by 
differentiating Equation 3-33 with respect to w and dividing 


by the mass m we obtain: 





ae ps 
Yes mages ——— (-C sine + CL COS o +G ws +@ sina) (3-34) 
Now, assuming “ small we express 

pl, S 

fw = Cu, +O) (3-35) 
where 
a: 

ante, e=6) 


nO) 





The quantity C,, represents the change in C, 
coefficient due to angle of attack changes. Graphically, it 


is well-known as "lift curve slope", as shown in Figure 3-3. 





Figure 3-3 


Lift Characteristic Curves for Different AR's 


The changes in “% , are linear functions of «, below 
the stall. Aspect ratio (AR), is defined simply as the ratio 


(span/chord) or in more general case as: 


(span)? b 
Ans cwing area ) 7 a (3-37) 


This ratio, determines the slope of the lift curve. As shown 
in Figure 3-3, low A%'s produce smaller lift coefficient G 
for the same angle of attack ~. 

The CG. derivative is always positive in the linear 
region. It is affected primarily by the wings, while the 
tail and fuselage contribute a little in aircrafts, but 


largely in missiles. Aeroelastic effects in addition may 
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cause bending or twisting of the wings that change the angle 
of attack and hence the value of C,. 
6. The Mw Derivative 

Is the rate of change in the pitching moment M, due 
to a w-velocity perturbation, divided by the Iy moment of 
inertia. 

oe al 
yy Ow (3-38) 

A w-velocity perturbation is associated with an angle 

of attack w which as we have seen is given by 
ays (3-39) 

When the angle of attack changes from its equilibrium 
position, the lift of the wing changes in the same sense, and 
that produces a positive or negative pitching moment 
depending on the center of gravity location. In addition, 
the lift of the tail changes and that contributes to the 
total pitching moment. These contributions together with the 
fuselage contribution make up this derivative. 

We desire a change in angle of attack to produce a 
pitching moment that tends to restore the original angle of 
attack. This requires a zero pitching moment in the 
undisturbed reference flight, a condition known as "trim 
COondge1ron.. 

The pitching moment, as we have seen, is expressed in 


terms of the non-dimentional pitching moment coefficient as 
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mi 2 


Because a small change in w-velocity doesn't affect 
the V velocity, by differentiating Equation 3-48 with respect 


to w we obtain: 





OM SOM oe 4 OM dye, 
Qw i  O« Ow vy Yo Wo P dol (3-41) 
and finally by letting V=U, we get 
_ pS bc 
Mw = 2 Tyy CHa (3-42) 


The center of gravity location plays a very important 
role in Cu, derivative. In fact Cu, is proportional to the 
distance between the center of gravity and the aerodynamic 
center, the center where the lift force acts. 

If the center of gravity coincides with the 
aerodynamic center, Cma is zero; if it is ahead of the 
aerodynamic center, Cu, is negative and the airframe is said 
to be statically stable; and if it is aft of the aerodynamic 
center, (Cm, is positive and the airframe is said to be 
Statically unstable. The more ahead the center of gravity is 
located, the steeper the Cm, curve becomes. All the above 
cases are shown in Figure 3-4. 

It is clear that a stable situation like the center 
of gravity at B or A, an increase in angle of attack from « 
to “ produces a negative pitching moment as Figure 3-4 shows, 


which tends to reduce the angle of attack. 
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Figure 3-4 


Cuy for Different Center of Gravity Locations 


Figure 3-4 also shows that for stable situations, for 
nearly small angles of attack, where the lift is nearly zero 
a mechanism must be provided to generate a positive “™. 
Here comes the effect of the tail which, being a negative 
angle of attack, generates a small lift force which though 
multiplied by a larger lever arm produces a balancing 
positive pitching moment. 

In other designs we find the tail being ahead of 


wings with positive angle of attack. This is a canard 


HOS 








configuration most usually found in missiles, also being the 
mechanism to generate a positive Cun. 

In Figure 3-4 Ce Curves are shown as linear 
functions of ow. This is not true for large variations in &, 
but can be nearly true for small variations around ™  . 

7. The Zy Derivative 

Is the rate of change in Z-force, due to a rate of 

change in w-velocity perturbation divided by the mass m 
7,-4.0 
y  ™ Ow (3-43) 


The relation &= w/Uo, leads us to the expression 


92 1 OZ 


ae) seeds: (3-44) 
According to Equation 3-16 the Z~-force is simplified 
to 
Z=W sd -L (3-45) 
Neither Wor 9 depend on rate of change in angle of 
attack o, since differentiating Equation 3-45 with respect to 


& we obtain 


ZA 
Ve DK (3-46) 
where from Equation 3-l 
OL 4 2 an 
age > 5, (3-47) 
Therefore, 
eee ok et te 8 
Lib = an Oy ale Be Tm De “mb 2 PY Fa (3-48) 
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To form a non-dimensional coefficient, we multiply 


and divide Equation 3-48 by c/2Uo so that letting V = Uo we 





~~ Men, Lik 
w ~~ md, 2 Pe? Wo DC uC/2Uo) (3-49) 
or _ pSe 
fw = - eta (3-59) 
where 9c 
eS 
Cla = BGae abe) (3-51) 


Ghiematanttty ch, represents a change in lift 
coefficient due to a variation in the rate of change of the 
angle of attack a. 

If the angle of attack changes rapidly, the pressure 
distribution on the wing, tail and fuselage, does not adjust 
itself instantaneously to its equilibrium value. This, a 
Change in lift and drag occurs to the sudden changes in angle 
Of “attack. 

Its sign is positive for low speeds but can be 
positive or negative on high speeds depending on aeroelastic 
effects. 

Also, a sudden change in angle of attack alters the 
downwash field as SEGA by the tail and a time delay is 
elapsed before the tail senses this alteration. 

All those effects are combined to give rise in the 


CH derivative. 
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8. The Mw Derivative 
Is the rate of change in the pitching moment M, due 
to a rate of change in w-velocity perturbation divided by the 


Jyy moment of inertia. 


M = pis oM 
"yy Ow (3-52) 


This derivative ®M/ow can be expressed as: 


OM 1 ‘OM 


—_ we ee 


Ww ~ Ue VK (3-53) 
According to Equation 3-3 
M = ie oV°Se a 
2 (3-54) 
Since V is independent of the rate of change in angle 
of attack «, by differentiating Equation 3-54 with respect to 


«x we obtain: 


‘aM pV*Sc ‘Xm 


om 2 Ve (3-55) 
Therefore, 
eee et ys, 
Tyy Ow  Velyy O* Uelyy 2 Oak (3-56) 


To form a non-dimensional coefficient, we multiply 


and divide Equation 3-56 by c/2Uo so that letting V = Uo, we 








get 9G. 
{ { 2 c 
My = a ples 7a if 
me! th «ICae/20,) oe 
Or — Sc? 
Mu = 4 Iyy Co (3-58) 





where 
9Cm 


“Mis * the vou) (3-59) 

Hor Quantity Cu y represents a change in pitching 
moment coefficient due to a variation in the rate of change 
of the angle of attack &. 

The change in lift that is produced by sudden changes 
in angle of attack discussed already, multiplied by 
corresponding lever arms, generates changes in pitching 
moment and that is why (Cm, arises. 

Same lag effects in the downwash field discussed, 
contribute to this derivative, as well as aeroelastic 
effects, that also disturb the air stream of the wing and 
tail, so that the pressure distribution over the aerodynamic 
Surfaces does not adjust instantaneously. 

The sign of this derivative is negative and is 
relatively important in longitudinal dynamics. Large 
absolute values of this derivative are desired because it 
cuuses damping. 

9. The Zq Derivative 

Is the rate of change in Z-force, due to a q-pitching 

velocity perturbation divided by the mass m 
ae U4 
4 om %q (3-68) 

According to Equation 3-16, the Z-force is simplified 

O's 


Z=Wos9 -L (3-61) 


Ibs: 





Neither W or G9 depend on the pitching velocity (q), 


Since differentiating Equation 3-61 with respect to q, we 


obtain: 
OZ OL 
(3-62) 


OL } 2 O. 
Therefore: 
‘Xr 
“9q (3-64) 


-£ 2 1 4 bh Lt tye 
Lat im 4  ™ % ” m Zev s 99, 


To form a non-dimensional coefficient we multiply and 


divide Equation 3-64 by c/2Uo so that letting V=Uo we get: 


Lh os —— 
La ~ Mm 2 p Us 5 20> OCq¢/2Ve (3-65) 
B)) & 
ss pSUec 
24 Aon Cog (3-66) 
where 
OC. 
(3-67) 


a 
i Cgc /2Vo) 

The quantity Chg represents a change in lift 

coefficient due to a variation in the pitching velocity. 


If the airframe pitches about its center of gravity, 


the angle of attack changes and a change in left of wing, 


tail and fuselage is developed that gives rise to this 


derivative. 
ale 





If has a small positive value in low speed flight but 
may have large values poSitive or negative in high speeds due 
to aeroelastic effects. 

Curved longitudinal flight paths give rise also to 
centrifugal forces on all components of the airframe that due 
to aeroelastic effects cause a bending of aerodynamic 
surfaces and, hence, a change in lift. 

In the past without computer solutions available, 
this derivative was considered zero to simplify the analysis. 

18. The Mq Derivative 

Is the rate of change in the pitching moment M, due 

to a pitching velocity perturbation, divided by the Iyy 


moment of inertia. 


Meat 2M 
+ Tyy 04, (3-68) 

According to Equation 3-3, 
M = = pV'Se Cy (3-69) 


Since V is independent of the pitching velocity q, by 


differentiating Equation 3-69 with respect to q, we obtain: 


OM 2 Cm 
2 *s, (3-78) 
Therefore, 
OM { 2 Ge 
My en, Cee ee 
4” Ty” ay P 2g (3-71) 


IL ALS 





To form a non-dimensional coefficient, we multiply 


and divide Equation 3-71 by c/2Uo, so that letting V=Uo we 


gets 
Reed U2 cM = 
Mq QIy, rs 2Uo OC 4c /2Vo) (3-72) 
OF 
_ pleS cz 
- = ~"s (sens) 
where 
Cc) 
$  OCqc /2ts) (3-74) 


The quantity (Mg represents a change in pitching 
moment coefficient due to a variation in the pitching 
velocity q. 

The change in lift that was produced, due to pitching 
velocity perturbations discussed already, multiplied by 
corresponding lever arms, generates changes in pitching 
moment and that is why Cg rises. 

The tail is the main contributor to this derivative 
which is also referred to as “pitching damping derivative". 
Aeroelastic effects and centrigugal forces generated due to 
curved longitudinal flight paths also contribute to this 
derivative. 

In longitudinal dynamics (Mg is very important. It 


has a negative value, and the larger it is, the more the 
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damping is increased; a condition that is highly desired. In 
high speed flights Cog may have positive or negative value 
depending on aeroelastic effects. It is also noted that 


tailless airplanes have a poor damping. 


3.3 LATERAL - DIRECTIONAL STABILITY DERIVATIVES 

Lateral-Directional stability derivatives are encountered 
in the lateral directional equations and, hence, they 
uniquely specify the Y-force, the L-rolling moment and the 
N-yawing moment. 

Dimensional analysis approach indicates that we can 
express the lateral-directional aerodynamic forces and 
moments in terms of their corresponding non-dimensional 


coefficients as: 


Y-force Y =C, 5 eV*S (3-75) 
L-rolling moment L=(C)3 pV’Sb (3-76) 
N-yawing moment N =C, pV*Sb (Ee TTe), 


where GY, (@% , & are the corresponding Y-force, rolling 
moment and yawing moment coefficients, and b is a reference 
length, usually taken the wing span. 
The Lateral-Directional Stability Derivatives can be 
evaluated as follows: 
1. The Yv Derivative 
Is the rate of change in Y-force due to a v-velocity 


perturbation divided by the mass m: 


a 


Way 





Y 


Y 


Vy = 7 (3-78) 


> 


A v-velocity perturbation causes a small change in 
the angle of sideslip B which is expressed according to 


Equation 2-133, approximately as: 


p ey es (3-79) 
so that 
ee OY 
1 = ae (3-80) 


Assuming that the V-velocity is independent of the 
angle of sideslip #, differentiating Equation 3-75 with 


respect to #, we obtain: 


ae 7 2°. > OB (3-81) 


Therefore, letting V=Uo we get: 


_ pos 
where 
Cy, = %4/ag (3-83) 


The quantity Cyp represents the change in the G side 
force coefficient due to a change in the angle of sideslip §. 
The main contributor of this derivative is the 
vertical tail. Secondarily the fuselage, and also the wings 


contribute a little. 
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Figure 3-5 shows the forces that the airframe 
experiences due to a v-velocity perturbation. First an angle 
of sideslip builds up and all aerodynamic surfaces are 


sensing a change in the approaching free stream. 
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Figure 3-5 


Forces Experienced Due to a v-Velocity Perturbation 


As a consequence, a force Fr due to the fuselage and 
a force Fyj due to the vertical tail appears. The wings also 
experience a different lift force. For the configuration 
shown in Figure 3-5, the right wing experiences more lift 
than the left wing because it senses a greater V-normal 
component of the approaching free Stream. Those lift forces 
may contribute to a small change in the Y-force. 

All these contributions make a rise to (, derivative 


which is very important in lateral-directional dynamics. It 





introduces a damping that tends to eliminate the angle Br 
because the resultant Y-force appears acting on a point aft 
of CG location and that generates the appropriate moment to 
do this. 

Its value is negative and it might be reasonable to 
desire a large value, for the # angle to die out quickly. 
Although high maneuverability airplanes must have a small 
negative value for this derivative, if they were to perform 
easy roll maneuvers and banked turns. 

2. The Ly Derivative 
Is the rate of change in L-rolling moment due to a 
v-velocity perturbation divided by the Ixx moment of inertia. 
> 
Tux Qv (3-84) 
It 1s also clearly expressed as: 
{ OL 
“Ulex 96 (3-85) 


Differentiating Equation 3-76, with respect to 4, we 


Ly 


obtain 
OL ! 2 aCe 
—=z=—pV°Sb —— 
3p 2 P 36 (3-86) 
Therefore, letting V=Uo we get: 
7 pU. Sb 
Ly = 9 Tex Ge (3-87) 
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where 
en 3 
@ ‘0p (3-88) 

The quantity Ge Eepresémres the change in the G] 
rolling moment coefficient due to a change in the angle of 
sideslip &£. 

It is usually referred to as the effective Dihedral 
Derivative. 

Dihedral is often used as a mechanism to improve 
lateral-directional stability. Figure 3-6(a) shows a head-on 
view of an airplane whose wings are tilted up at some 
dihedral angle with respect to the horizontal. 

We have already discussed that due to a v-velocity 
Doacuroatnaon, according to the configuration shown in 
Figure 3-5, where a positive sideslip angle £ is shown, right 
wing generates more lift than the left. This is true for the 
Swept wings. If the wings are not swept, the mechanism to 
produce the same effect is the dihedral. Usually even swept 
wings are given a Small dihedral to improve 
lateral-directional stability. 

To understand the dihedral effect, consider the 
aircraft shown in Figure 3-7 with non-swept wings but with 
some dihedral angle. The aircraft is distributed by an angle 


of sideslip A. 
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Figure 3-6 


Dihedral Effect on Lateral- SW BNE Oe Stability 
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Figure 3-7 


Dihedral Effect on Non-Swept Wings 
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The perpendicular components of the free stream 
velocity are equal but the angle of attack that the right 
Wing projects to this perpendicular component is greater than 
the angle of attack that the left wing projects, as an effect 
of the dihedral angle. 

As a consequence, the right wing generates more lift 
than the left wing aS shown in Figure 3-6(c). Those lift 
forces contribute to a change in the rolling moment 
coefficient. 

Figure 3-6(d) and (c) show the effect of the 
high-wing placement which has the same effect as the dihedral 
angle. It is clearly shown that in the low wing placement 
configuration, the fuselage does not permit the vertical 
component of the free stream to act upon the right wing to 
produce a stabilizing moment, whereas in the high wing 
placement configuration the fuselage effect is nearly zero 
and that makes a stabilizing rolling moment to show up. 

In addition, the force acting on the vertical tail 
Shown in Figure 3-5, multiplied by the corresponding vertical 
arm relative to the x-longitudinal axis, contributes also to 
the same effect. 

Pg is nearly always negative in sign and the 
corresponding rolling moment L, is negative for a positive 
sideslip angle f. It is a very important derivative in 


lateral directional dynamics because it causes a damping and 
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improves the stability. Although high negative values of 
this derivative will tend to decrease the maneuverability of 
the airplane in roll movements. In missiles a low negative 
value of this derivative is desired. 
3. The Nv Derivative 
Is the rate of change in N-yawing moment due to a 


v-velocity perturbation divided by the Izz moment of inertia. 


Tez Ov (3-89) 
It is also clearly expressed as; 


LN 
tee aye (3-90) 


N 
Differentiating Equation 3-77, with respect to $, we 
obtains: 
ON 2 OCn 
ye 2h Se oe (3-91) 


Therefore, letting V=Uo we get: 


_ pU.Sb 
al Vin 6% (3-92) 
where 
Cng = Cn /Og (3-93) 


The quantity (mg represents the change in the G 
yawing moment coefficient due to a change in the angle of 


Sideslip 8. 
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It is usually referred to as the Weathercock 
Derivative. The main contributor of this derivative is the 
vertical tail and secondarily, the fuselage. If the airframe 
experiences a v-~-velocity perturbation and an angle of 
sideslip & appears, the free relative stream strikes the 
Fuselage and the vertical tail in such a manner that the 
force Fr , shown in Figure 3-5, produces a destabilizing 
yawing moment and the force fvt a stabilizing moment. Wings 
contribute a little in the positive (stabilizing) sense. The 
total value of this derivative is positive, signifying static 
directional stability. 

The (mn, derivative, is very important in 
lateral-directional dynamics. The vertical tail must be 
properly designed to provide a sufficient weathercock effect. 
At large, sideslip angles a stall may occur due to flow 
separation and a catastrophic sideslip divergence may result. 
A low aspect ratio prevents stalling but decreases stability. 

Usually a large positive value of the Ong derivative 
is desired, by extending the vertical tail. A large value 
helps the pilot in good coordinated turns and prevents 
Sideslip. 

4. The Yy Derivative 
Is the rate of change in Y-force, due to a rate of 


change in v-velocity perturbation, divided by the mass m. 


a? 


a 


Yo = Aye (3=94) 


3|> 


eZ 





It is also clearly expressed as: 


7, aes 
" "Vom 98 (3-95) 





Differentiating Equation 3-75 with respect to B we 


obtain: 


OY LL ytg 2& 
me 2 PY S % (3-96) 


To form a non-dimensional coefficient we multiply and 


divide Equation 3-96 by b/2Uo. 


oY Lt pytg &. _% 
me 2° > au IA) (3-97) 


Therefore, letting V=Uo we get: 


_ pSb 
Yo am S46 (3-98) 
where 
_ oCy 
Ye 3Cu 7am (3-99) 


The quantity Cy, represents the change in the 
coefficient with variation in rate of change of the sideslip 
angle 8. 

If the sideslip angle # changes rapidly, the pressure 
distribution on the vertical tail and fuselage does not 
adjust itself instantaneously to its equilibrium value. 

This effect, as well as aeroelastic effects, give 
rise in Cys derivative but its value is very small and is 


usually neglected. 
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5. The Ly Derivative 
Is the rate of change in L-rolling moment, due to a 
rate of change in v-velocity perturbation divided by the Ixx 


moment of inertia. 


La = de . OL 
alex eal’ (3-108) 





It is also clearly expressed as: 


eee le 
eo (3-101) 





Differentiating Equation 3-76, with respect to a, we 


obtain: 


oC 


OL 7 f 2 4 
c= v Soar (3-162) 
To form a non-dimensional coefficient 
OL 4 2 b pa 
ee eC Se 
of 2 P ? 2 Uo OC Bb /2Ve) (3-193) 


Therefore, letting V=Uo 





 psh 
-* 4 Ixx Ce i (3-184) 
where 
Ce ‘Ce 
Ci OC Bb/ 2) (3-105) 


The quantity Ce represents the change in the G 
rolling moment coefficient with variation in rate of change 


of the sideslip angle &8. 
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The change in the side forces on the vertical tail 
and the fuselage produced due to sudden changes in the 
sideslip angle #, discussed already, multiplied by 
corresponding lever arms provide a total change in the Co 
rolling moment coefficient. 

This derivative has a very small value because the 
lever arms are very small and it is usually neglected. 

6. The Ny Derivative 

Is the rate of change in the N-yawing moment, due to 
a rate of change in v-velocity perturbation divided by the 
Izz moment of inertia. 

Wet Ss : 
le ov (3-196) 


It is also clearly expressed as: 


ee 
i iss oye (3) 
Differentiating Equation 3-77, with respect to By we 
obtain: 
ON { 2 oCn 
a = — Sb 7 ae ae = 


To form a non-dimensional coefficient: 


ON 2 


Therefore, letting V=Uo we get: 


b Wy 
QU. DC fb/2U.) (3-199) 





Ny = ae Cn (3-118) 





where 
C5, a eae 
6 9(Bb /2U.) (3-111) 

The quantity (ng represents the change in the G 
yawing moment coefficient with variation in rate of change of 
the sideslip angle §. 

The change in the side forces on the vertical tail 
and the fuselage produced due to sudden changes in the 
sideslip angle &, discussed already, multiplied by 
corresponding lever arms provide a total change in the Cn 
yawing moment coefficient. 

The lever arm of the vertical tail this time is long 
and hence, this derivative may be considerable. Positive 
values of this derivative increase the damping in the yawing 
motion and are desirable in general. 

7. The Y, Derivative 

Is the rate of change in the Y-force, due to a 
p-rolling velocity perturbation divided by the mass m. 

A 
nn (3-112) 
Differentiating Equation 3-75, with respect to p, we 


obtain: 


(3-113) 


Wee 





To form a non-dimensional coefficient 


“OY 


Nt ayee boy 
Op 2? > ay Tpb/2v) (3-114) 


Therefore, letting V=Uo 


p Uo Sb 
— Am CYp (3-115) 
where 
6, - 2% 
Yp = OCpb/2vp ) (3-116) 


The quantity a. represents the change in the 
Y-force coefficient with variation in the rolling velocity p. 

When an aircraft has a roll velocity p, a force is 
acting on the vertical tail that usually opposes the roll 
movement. This consists the main contribution to this 
derivative and is small enough to be neglected. 

The force acting on the vertical fin may be affected 
by the sidewash from the wing that gives a slightly different 
direction of the airstream. 

8. The L, Derivative 

Is the rate of change in the L~-rolling moment, due to 

a rolling velocity p perturbation divided by the Ixx moment 


of inertia. 


OL 
Lp ace OD (3-117) 
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Differentiating Equation 3-76, with respect to p, we 


obtain: 
ao pV?Sh 
Op 2 ‘Op (3-118) 
To form a non-dimensional coefficient, 
Se 
a) Z 20 (pb /2Uo) (3-119) 
Therefore, letting V=Uo 
) p Up SY C 
fe ins - Pp (3-128) 
where 
Ce = Xe 


P —  QCpb/2U.) (3-121) 

The quantity 4, represents the change in @ rolling 
moment coefficient with variation in the rolling velocity p. 

This derivative is also referred to as "roll damping 
derivative", because it has usually a negative or damping 
effect coming mainly from the wings and secondarily, from the 
vertical and horizontal tails. 

When the aircraft rolls, the force acting on the 
vertical tail, shown in Figure 3-8, multiplied by the 
corresponding lever arm, produces an opposing rolling moment. 
In addition, on the down-going wing, the effective angle of 
attack is increased and on the up-going wing, decreased. 
This makes the lift on the down-going wing greater and thus 
an opposing rolling moment is produced. Similar effect, on a 


smaller scale, is produced by the horizontal tails. 
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down - qving wing up- going wing 


Figure 3-8 


Forces Due to Rolling Velocity (p) 


This is true as long as the down-going wing is not 
Flying near the stall due to a large angle of attack 
experienced in some cases. In such caseS a positive effect 
will result in autorotation. 

The derivative Ce, Wommaruite important @en 
lateral-directional dynamics. The value of this derivative 
directly affects the ailerons design. A small negative value 


of this derivative will result in a better aileron input 
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response while a large negative value will make the airframe 
heavily damped in roll maneuvers. 
9. The No Derivative 
Is the rate of change in the N-yawing moment, due to 
a rolling velocity p perturbation, divided by the Izz moment 


of inertia. 


{ ON 


——=——_ e en 


P Tee % (=i) 
Differentiating Equation 3-77, with respect to p, we 
obtain: 


ON 1 2 ‘OGn 
Re 2 es’ Op Belge) 





To form a non-dimensional coefficient, 


ON 2 b ner 
eee | ee ee 
Op 2 P a 20, O(pb/2B) (3-124) 


Therefore, letting V=Uo 





_ pUSb’ 
Np = Tee Cn (32125) 
where 
C = BGs 
™e  OCpb /2be) (3-126) 


The quantity ("p represents the change in Gn yawing 
moment coefficient with variation in the folling velocity p. 

This derivative usually has a negative or positive 
effect coming from the wings, the vertical tail and 


secondarily, from the horizontal tail. 
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When the aircraft rolls, the force acting on the 
vertical tail, shown in Figure 3-8, multiplied by the 
corresponding lever arm produces a yawing moment towards the 
same side of the roll motion, i.e., roll to starboard causes 
yaw to starboard. This contributes a positive effect and 
depands upon the vertical tail geometry and the Ssidewash from 
the wings. 

In addition, the contribution of the wings can be 
considered to be made of two parts, one opposite the other, 
the total effect being negative. When the airframe rolls, on 
the down-going wing, the effective angle of attack increases 
and on the up-going wing decreases, as shown in Figure 3-8. 
As a consequence, the lift and drag are greater on the 
down-going wing than on the up-going wing. The projections 
of the drag forces on the xoy plane usually produce a yawing 
moment towards the same side of the roll motion. This 
produces a positive effect while the projections of the lift 
forces produce a yawing moment towards the opposite side of 
the roll motion, and this produces a negative or damping 
effect. 

The overall value of this derivative is usually 
negative, i.e. roll to starboard produces a yaw to port. At 
higher speeds the numerical value decreases and may become 


positive. 
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Tas ewe riveatweve is fairly important in 
lateral-directional dynamics, as far as roll maneuvers and 
their effect to yaw. 

16. The Yr Derivative 

Is the rate of change in the Y-side force, due toa 

yawing velocity perturbation r, divided by the mass m. 
| 
Game Ga) 


Differentiating Equation 3-75, with respect to r, we 


obtain: 


‘OY ! 


v7 2 Or (2223) 


To form a non-dimensional coefficient, 


eoers 2S 
7 ~2P" > Dw Beb/2D) (3-129) 


Therefore, letting V=Uo 


pUo$ b 
pe oe Cie (3-138) 
where 
eeaeneeees 
Ur = DCrb yay (3-131) 


The quantity ‘4, represents the change in Cy Y-force 
coefficient with variation in the yawing velocity r. 
When the aircraft has a yaw velocity r, as Figure 3-9 


shows, a force is acting on the fuselage that usually has a 
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negative or damping effect on the yaw motion and a force is 


acting on the vertical tail that usually has the same effect. 





Figure 3-9 


Forces Due to Yaw Velocity (r) 


Both contributions are small enough to be neglected, 
therefore, this derivative is of little importance in 
lateral-directional dynamics. 

11. The §, Derivative 
Is the rate of change in the L~rolling moment, due to 


a yawing velocity r perturbation, divided by the Ixx moment 


of inertia. 


ye CISA 
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Differentiating Equation 3-76, with respect to r, we 


obtain: 


OL J 2 ol | 
foe PY Sb Bash 


To form a non-dimensional coefficient, 


OL A 2 b ‘Ce 
or a7 pv Sb We | 9Crb/2Vo) (3-134) 


Therefore, letting V=Uo 


p Ub S b 
ee — e 
4 Ixx r (3-135) 
where 
G - OC» 
Yr - Orb /2Uo) (3-136) 


The quantity @, represents the change in Ce rolling 
moment coefficient with a variation in the yawing velocity r. 

When the aircraft is yawing about its vertical axis, 
the outer wing moves faster than the inner. As a 
consequence, the lift on the outer wing is greater than that 
of the inner and, hence, a rolling moment is produced towards 
the same side of the yaw motion. This contributes a positive 
effect, i.e. yaw to starboard produces roll to starboard. 

In addition to this major wing contribution, the 
Force acting on the vertical tail multiplied by the 


corresponding lever arm, produces the same rolling effect. 


Le 








The total value of this derivative is usually 
positive, but it iS not so important in lateral-directional 
dynamics. 

12. The Nr derivative 

Is the rate of change in the N-yawing moment, due to 
a yawing velocity r perturbation, divided by the Izz moment 
of inertia. 

\. Sp BL 
Izz Or (3-137) 


Differentiating Equation 3-77, with respect to r, we 


obtains: 
ON A 2 ‘oCn 
me 7°" >? > (3-138) 


To form a non-dimensional coefficient, 


ON { 2 b On 
Gp 2 did 20, rb /2Uo) (3-139) 


Therefore, letting V=Uo 


2 
sone PS Be 


Mig tt (3-149) 
where 
OC 
Cn, = = 3-141 
: O(rb/2Zlb ) ( 41) 


The quantity (n, represents the change in CG yawing 
moment coefficient, with a variation in the yawing 


velocity rc. 
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This derivative is known as the "yaw damping 
derivative", because it usually has a negative effect. 

The main contributor is the vertical tail. When the 
aircraft 1S yawing, the force acting on the vertical tail 
multiplied by the corresponding long lever arm, produces a 
considerable yawing moment that opposes the initial yawing 
disturbance. 

In addition, the drag on the outer wing is greater 
than that of the inner wing. That produces a yawing moment 
that opposes the initial yawing disturbance. The force on 
the fuselage also contributes an opposing effect. 

Miiesmeaerivative is very important in 
lateral-directional dynamics. Its total effect is negative 


and large values are desired for more effective damping. 


3.4 CONTROL SURFACE DEFLECTIONS 

We have discussed so far all the important stability 
derivatives that rise due to the airframe geometry. They 
uniquely specify the forces and moments about the center of 
gravity provided that they are multiplied by the correct 
Factor that depends upon the airstream and airframe 
characteristics. 

A stability analysis is possible now for determining the 
longitudinal and lateral directional modes of motion, 
corresponding to an initial disturbance, i.e. homogeneous 


solution to the differential equations. 
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In addition, some means of control in longitudinal and 


lateral directional motion are available to the human or 


automatic pilot. They are shown in Figure 3-19. 


os elevator control 7 





rudder coutrol 


om @& 


Figure 3~10 


Basic Longitudinal and Lateral Control Systems 





Elevator, controls the longitudinal motion of the 
airframe providing nose up or nose down moments as a response 
to corresponding control surface deflections. 

Aileron and rudder, control the lateral directional 
motion providing rolling and yawing moments correspondingly 
as a response to the corresponding control surfaces 
deflections. 

The symbols used for the corresponding angles of rotation 


and control surface deflections are shown below: 


[a 
& 6 W 
p . 
angle of rotation ( pitch) ( yaw) 
5 mn 5 
oa aia a ah 


In the longitudinal and lateral directional equations of 






motion, Equations 2-152 and 2-153, respectively, we have 
denoted by the symbol 6 the control input due to control 
Surfaces deflections. 

The right-hand side of the longitudinal Equation 2-152 is 


expanded as follows: 


X§ = Xm 
Z5 = Zn (3-142) 
Ms = Mn+Mn 
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and the right-hand side of the lateral directional 
Equation 2-153 as; 
VS Bais 
Lae Ee Sees (3-143) 
Nj = Ne +No+Ne 

In addition, to the right-hand side terms that take care 
of the main control surfaces, there may exist terms rising 
from the auxiliary surfaces such as flaps, stabilizers, dive 
brakes, etc. 

A general approach suggests to examine the subscripted 
derivatives instead of the particular derivatives dedicated 
for a specific control surface deflection, i.e. the total 
effect by all the possible combinations of the control 
surface deflections. The approach is Similar to that 
followed in Sections 3.2 and 3.3. 

1. The %§ Derivative 

Is the rate of change in the x-force due to a 


control surface perturbation divided by the mass m. 


re) 
S| 


X53 (Jaan 


= 
am 
By 6 we denote a general control surface deflection. 
For longitudinal control it might be elevator, stabilizer, 
flap, slat, dive brakes, etc. 
According to Equation 3-27, the X-force is simplified 
EO: 


{ 2 
X=-57 pV SG +T (3-145) 
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We can assume that neither V or T depends on the 
control deflection. Thus differentiating Equation 3-145, 


with respect to , we obtain: 


ny ae OS (3-146) 
If we let V=Uo and 


Co, = 96/98 (Sia) 


OX I /* S IG 


we arrive at the expression: 


p Vo" S 
Die CDs (3-148) 


Xs = - 

The quantity Qs represents a change in the drag 
coefficient due to a variation in a longitudinal control 
surface deflection. 

A control deflection, positive or negative, increases 
the drag due to the greater projected frontal area, but this 
sometimes is desirable because it provides damping in an 
oscillatory perturbed motion. 

2. The 23 Derivative 


Is the rate of change in the Z-force due to a } 


control surface perturbation divided by the mass m. 


£5 im 3 (3-149) 
According to Equation 3-16, the Z-force is simplified 


COs 
Zea Wicos Ore |e (3-156) 
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or 
{ 2 

= - — Se 

Neither W or 9 depends on the d control deflection. 


Thus differentiating Equation 3-151, with respect to d we 


obtain: 
CA eee 2G. 
ae ES Benao) 
If we let V=Uo and 
Cus = C, / 9S (3-153) 
we arrive at the expression 
2 
_ ps 
A ocean Oe (3-15 4) 


The quantity Ci represents a change in the lift 
coefficient due to a variation in a longitudinal control 
Surface deflection. 

This derivative is a measure of the effectiveness of 
the elevator mainly, in producing a change in the lift force 
acting on the tailplane. A positive elevator deflection 
produces an increase of lift, i.e. negative Z-force. 

Tt has a very small value on aircrafts with tail but 
a large value on tailless aircrafts. 

3. The MS Derivative 

Is the rate of change in the M-pitching moment due to 
a control surface perturbation divided by the Iyy moment of 
inertia 


| °M 
eee (3-155) 
M9 dyy res) 


144 





According to Equation 3-3, the pitching moment is 


expressed as; ; 


{ 
M= 5 pVScCy (3-156) 
Differentiating Equation 3-156, with respect to 5 we 
obtain: 
oM ae pV*Sc Xm 
90 82 0d (3-157) 
If we let V=Uo and 
Cus 2 Om /9g (3-158) 
we arrive at the expression 
acre 
Me 2 pe Se 
8 QIyy CM § (3-159) 


The quantity Cus represents a change in the pitching 
moment coefficient due to a variation in a longitudinal 
control surface deflection. 

It is commonly referred to as the elevator power 
derivative because it is a meaSure of the elevator control 
effectiveness. 

The change in Z-force produced by the elevator 
control surface deflection, multiplied by a rather long arm, 
provides a pitching moment in the negative sense, i.e. a 
positive elevator deflection results in a nose-down 
(negative) pitching moment. 

That is the primary function of the elevator and is 


considered to be the most important of all the control 
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surface functions. It adjusts the angle of attack and it 
trims the vehicle in its longitudinal motion, by providing 
Opposite pitching moments to accommodate disturbing moments 
and center of gravity relocations. 

For a design point of view, the larger the center of 
gravity range is, the greater the elevator pitching moments 
must be provided, i.e. larger (ms values. 

4. The YS Derivative 
Is the change in the Y-side force due to 2 d control 


surface perturbation divided by the mass, m. 


q@ 


aL 


| 
Y5 = m 95. (3-168) 


Differentiating Equation 3-75, with respect to é we 


obtain: 
ote V2s oy 
3 2? 25 (3-161) 
If we let V=Uo and 
ys = Gy /25 (3-162) 
we arrive at the expression: 
Ye = A> Cy; (3-163) 


The quantity ‘y; represents the rate of change in the 
Y-side force coefficient due to a variation in a lateral 


control surface deflection. 
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This derivative rises mainly from the rudder. A 
positive rudder deflection is defined to produce a positive 
side force. Therefore, the value of this derivative is 
positive but small enough to be neglected. 

Rudder deflections also cause an increase of the 
total drag of the vehicle but that has a negligible effect on 
the longitudinal motion, and is neglected. 

Ee etne ls Derivative 

Is the rate of change in the L-rolling moment due to 

a ) control surface perturbation, divided by the Ixx moment 


of inertia. a 


xx | 98 (3-164) 


Differentiating Equation 3-76, with respect to d, we 


obtain: 

OL { 2 ‘Xp 

sesh’ S bans (3-165) 

If we let V=Uo and 
Ce, = IG /o5 (3-166) 

we arrive at the expression: 

a OUs Sb ; 

Manor, (Gaalogn 


The quantity (; represents the change in the © 
rolling moment coefficient due to a variation in a lateral 


control surface deflection. 
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The side force generated by a rudder deflection, 
produces a rolling moment in the positive sense, i.e. a 
positive rudder deflection generates a positive rolling 
moment. In airplanes, with large vertical tails, this effect 
is negligible and is often neglected. 

The main contribution to this derivative is from the 
ailerons, and it iS a meaSure of the aileron effectiveness. 
A positive @€lleron deflection is defined as the one to 
produce a positive rolling moment, i.e. left aileron down and 
right aileron up, therefore, the value of this derivative is 
positive. 

In lateral-directional dynamics, this derivative is 
probably the most important because it establishes the 
maximum available rate of roll of an airplane which is a very 
important consideration in combat environment. 

It also plays an important role in take-offs, and 
landings, where adequate and rapid lateral control is highly 
desirable. 

Rapid rate of change of aileron deflections gives 
also a rise to this derivative due to lag effects, 1.e. not 
immediate adjustment of the airstream. 

6. The Nj Derivative 

Is the rate of change in the N-yawing moment due to a 

3 control surface perturbation, divided by the Izz moment of 


inertia, 
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cx oe ac 
Tas (3-168) 


Differentiating Equation 3-77, with respect too, we 





ebtain: 
ON | 2 Qe 
ame TT PY SS a | (3-169) 
If we let V=Uo 
Cus = Wn/ 95 (3-178) 
we arrive at the expression: 
N : pUs Sc ‘é 
oar SDS bees (3—171) 


The quantity n, represents the change in the G 
yawing moment coefficient due to a variation in a lateral 
control surface deflection. 

The side force generated by a rudder deflection, 
produces a yawing moment in the negative sense, i.e. a 
positive or to the left rudder deflection, produces a 
negative yawing moment. This derivative is a measure of the 
rudder effectiveness. 

Rapid rate of change of rudder deflections gives also 
a rise to this derivative due to lag effects. 

This derivative plays an important role in crosswind 
takeoffs, and landings, and counteracts adverse yaw in 


rolling maneuvers. 
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Finally, aileron deflections give also a rise to this 
derivative since each aileron introduces a drag force of 
different magnitude. This effect may be positive or negative 
depanding on the rigging of the aileron and the angle of 


attack of the airframe. 


3.5 DETERMINATION OF THE STABILITY DERIVATIVES AND OTHER 
CONSIDERATIONS 

All the stability derivatives discussed so far rise from 
the flight vehicle's geometry, being highly dependent upon 
the design layout. We have referred mostly in this chapter 
to the conventional airplane with horizontal tail and we have 
seen that the derivatives are made up of a number of separate 
contributions. 

We had strictly remained in the dimensional form of the 
Stability derivatives and we simply expressed them in terms 
of airstream characteristics, airframe geometry and 
non-dimensional derivatives. 

Dimensional form has the advantage of obtaining real time 
solutions rather than working with normalized time mass and 
inertia components that must also be non-dimensional in 
considering that form. 

The advantage of non-dimensional derivatives, i.e. 
independent of size and airstream characteristics 1s 


maintained because dimensional derivatives are expressed with 
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the corresponding non-dimensionals, some of which though do 
depend upon the airstream and airframe characteristics. 

Thus, some of the dimensional expressions and some of the 
non-dimensionals are constants, but they can be expressed in 
polynomial form. 

As an example, for the METEOR 7: 

G, = -0@.052 - 9.999 ©, 
Co. = 6.0295 + 9.159 ¢, 
(np = -§9.942 C, 

(n, = -9.182 - 9.052 C, 

For others, plots may exist indicating their variation 
with respect to © or mach number. 

During the development of an aircraft, it is important to 
Know the numerical values of the derivatives under various 
conditions of flight. In the first stage, theoretical 
estimates can be made supported by wind tunnel tests and 
model experiments. 

By then, Simulation in digital computers can be made to 
obtain response characteristics to initial disturbances or 
control deflections and determine whether the values of the 
derivatives are satisfactory. 

When the aircraft has been built, flight tests are 


carried out to confirm the theoretical with the actual 


responses. : 
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There also exist flight test techniques for determining 
the values of the derivatives. 

Tables 3.1 and 3.2 summarize the longitudinal and lateral 
directional derivatives including definition, unit, 


expression and auxilliary relations. 
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TABLE 3-1 


LONGITUDINAL STABILITY DERIVATIVES 
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TABLE 3-2 
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LATERAL DIRECTIONAL STABILITY DERIVATIVES 


Quantity | Definition 


Aux. Relation 
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CHAPTER 4 


LONGITUDINAL DYNAMICS 


4.1 INTRODUCTION 

Having specified all the vehicle's dynamic properties, we 
will now examine the dynamic motion, i.e. the vehicle's 
response to disturbances and control inputs. 

Since any flight vehicle is a dynamic system, described 
by differential equations of time, we must expect in a 
particular situation, oscillatory or exponential converging 
or diverging response. 

Any particular response will be associated by its 
corresponding characteristic root location in the complex 
plane, as we have clearly seen in the dynamic response of the 
Spring mass damper system. 

The characteristic polynomials in the Laplace variable 
(s) involved, are usually of the third degree or higher and 
thus hard to solve by hand. 

In the past, a lot of approximations were introduced and 
even reduction of airframe's degrees of freedom technique was 
required to analyze coarsely a particular mode of motion, by 
reducing in effect the degree of the characteristic 
polynomial. 

In this chapter, we will analyze the longitudinal motion 


of a conventional airplane. We will first discuss the 
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Laplace Transform solution by obtaining the roots of the 
associated characteristic polynomial. Then we will apply the 
reduction of degrees of freedom technique to isolate and 
analyze a particular mode of motion. Finally, the equations 
of motion will be rearranged in a form that is more suitable 
for computer solutions. 

As an example, the longitudinal dynamic motion of the 
aircraft shown in Figure 4-1 will be analyzed with the 


parameters indicated. 


4.2 CHARACTERISTIC POLYNOMIAL OF LONGITUDINAL MOTION 
Longitudinal equations of motion are farther simplified 
by neglecting the effect of the Zq and Zw derivatives. 


Experimental results have shown that the effect of neglecting 


them is only of secondary importance. Further a straight 
level flight is assumed as the reference flight condition 


(1.e. ¥o =9), hence, Equation 2-152 can be reduced to: 


-$-Xy -Xw O g v X§ 

-Zy -Zw -U, oy; wt] 48 5 

-My  -sMw-Mw s-Mq O 1 E M5 (4-1) 
0 0 { set, Gcieeeo 


The characteristic polynomial is found by evaluation of 
the above determinant.: The Equation 4-1 can be reduced in 3 
equations of the original %,w,q, longitudinal variables by 


avoiding to write down the definition equation for q, but 
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altitude (4) 20000 
weight Clb) 30500 
mach number 0.638 
airspeed CAlsec) G60 





X%y =--0.0097 ty < ~0,0955 Uy a 0.0 
Xy = 0.0016 2w= - 1.43 Mw = -0.0235 
Xs = 00 25= -C8 Uw = -0,0013 
Ma ~ - 1,92 
U3 = = 26.10 
— 53.375 a 56.167 
17.583 
<2) 
(all dimensions in ft) 
Figure 4-1] 


Conventional Airplane and Parameters 
Used for the Numerical Example 
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that doesn't simplify the hand calculations. The determinant 


comes out to be: 


: fees es Gs 4 Ds +E =O (4-2) 
where 

al! (4-3) 

B= = Mg -UoM ws - Zw -Xy (4-4) 

C= ZyMq -UoMw -XwZv + Xv (Mg +UoMw +Zw) (4-5) 

D = -Xv (Zw Mq -UoMw) + Zu (Xw Mg +g Mw) - Mu (UoXw-g) (4-6) 

E = 9 (Zo Mw - Mu Zw) (4-7) 


Equation 4-2 is referred to as the longitudinal stability 
quartic. The roots in this case for any conventional vehicle 
come out to be complex with negative real parts, which are 
associated with two oscillatory convergent motions. 

The overall longitudinal motion of the vehicle is a 
Superposition or a combination of these two oscillatory 
modes. 

One of the modes usually is associated with a pair of 
complex roots located near the origin of the negative complex 
plane. This mode is called "Phugoid motion" and is a fairly 
lightly damped oscillation with long periodic time. 

The other mode uSually is associated with a pair of 
complex roots located far from the origin of the negative 
complex plane. This mode is called "Short period motion" and 
isa heavily damped oscillation with very short periodic 


time. 
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The longitudinal oscillatory modes are 


Figure 4-2. 


short pertod oscillation 
/ far oscillation 
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Figure 4-2 


The Longitudinal Oscillatory Modes 


For the given example we find 


A= 1.8 

B = 4.2177 

C = 18.2962 
D= 8.1814 

E = 8.0722 


and the characteristic polynomial is 
544 42177 s? + 18.2962 57+ O1814 5 + 0.0722 =O 
which has roots 
-8.0045 + j 8.9627 
—-2.1043 + j 3.7184 


is), 


shown in 


(4-8) 
(4-9) 
(4-18) 
(4-11) 
(4-12) 


(4-13) 


(4-14) 
(4-15) 





The pair of roots -9.80945 + j 9.8627 correspond to the 
phugoid motion and the pair of roots -2.1043 + j 3.7184 
correspond to the short motion. 

The fourth degree polynomial can be factored in two 
second degree polynomials. Each one of the second degree 
polynomials is associated with the corresponding mode. 

Further it is more convenient to put each one of the 
polynomials in the familiar form 

s?+ 25unS + Wy (4-16) 
so that the natural frequency W,, the damping 5 are figured 
out by inspection while the damped frequency can be 
determined by: 


1/ 
Wy = Wn (4-5?) : 


(4-17) 
If we denote with the subscript (SP) the short period 
quantities and with the subscript (PH) the phugoid 
quantities, the characteristic polynomial may be written in 
the form: 
(s?+ 25,, Wn, s +Wn., )(s?+ 25,, wn,, 5 tWn,,) =O (4-18) 
Those quantities and their relation is shown in the 
complex plane in Figure 4-3. 
The roots of the characteristic polynomial are clearly 
Shown in Figure 4-3 to be 
= S59 ieee: § Lox (4-19) 


~ Ose Wnsp = j Ase ea) 
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Or 


., 2 
“Seu Way + j Wnpy (1- Spy) | (4-21) 


2 


- Sse Wisp + | Wn sp (1-5 sp) 


(4-22) 


oY 
Wd, sWn,, TS Me 






Wal a, e WM oi, G-Sa yy 


f - Sse Wns» 
- Sph Yn ph 


~——-- oe af 


Mwerrrerr rrr 


Figure 4-3 


Longitudinal Roots and Characteristic Quantities 
in the Complex Plane 


Other quantities of importance in the longitudinal 
dynamics are: The phugoid period denoted by Ty, , the short 
period period T., given by: 

Tey = 2a /wy,, Csec) (4-23) 
ss 3 AD ieee Csec ) (4—24) 
and the time for half amplitude for phugoid and short period 


given by: 


fhe) 3h 





0.693] 


tM pu Goa Wn (4-25) 
4 0.6934 
"esp ~ Ssp Wnse ame) 


The Equations 4-25 and 4-26 come from the fact that the 
oscillatory response decays exponentially described by the 
envelope pee where Re is the real part of the 
characteristic root. Therefore, to obtain half of the 
amplitude the relation 

ae (4-27) 
must exist. By taking the natural logarithm of both sides we 
obtain the above relations. 

For our example, the characteristic polynomial of 
Equation 4-13 can be factored in the form 


(s* + 6.0999 s + 0.090395) (s? + 4.2086s + 18.2546) (4-28) 


From this, by inspection we can determine: 








Uno, = (6.90395) * = 6.0628 rad/sec (4-29) 
Wy = (18.2546)"2 = 4.2725 rad/sec (4-38) 

0.0090 

O717 _ 

Spa 9x 0.0628 y ieee) 

4.2086 
Sees eeaeres (°7925 eee! 

and also 
0.6934 

Ly oy ~ OonI7 + 0.0628 — os (4-33) 
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0.694 
eet = (OSS) wee 


ee ee 

"esp ~ 0.4925 « 4.2725 (4-34) 

Ou,, = 0.0628(1-6.07172)"2 = 9.9627 rad/sec (4-35) 

Wisp = 4.2725(1-6.49257)* = 3.7184 rad/sec (4-36) 
an 

Tew = Doe? 400.2! sec (4-37) 

7 2n 

ise = 37gq = 169 sec VEY, 


It 1S worthwhile to note that: 
(a) The damped frequency is indicated by the 
imaginary part of the corresponding characteristic root. 
(b) The real part of the characteristic root 
indicates the exponentially decaying envelope of the 
corresponding oscillatory motion, i.e. the damping. 


(c) The following relations hold 


B = 254 Wn,, + 25sp Wn sp (4-39) 
C= Wp, +45 py Wry dsp Ynep + Wee (4-49) 
D = 254 Wnpy Wrp + 25sp Wn, WA oy (4-41) 
E = Wipy Wigs (4-42) 


The two types of dynamical longitudinal oscillations are 
Shown in Figure 4-4, 

The forward velocity v and the pitch angle 9 are usually 
excited in the phugoid mode while the w velocity and, hence, 


the angle of attack % remains nearly unchanged. 
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Maxime Seepol 
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Short- period lonqitudmeai oscillation 


Figure 4-4 


Dynamic Longitudinal Oscillations 


In contrast, the angle of attack is heavily excited in 
the short period mode as well as the pitch angle 9, while the 
forward velocity remains nearly unchanged in that mode. 

The pilot can generally control the phugoid oscillation 
although because of his slow reaction it is possible to 
worsen to unstable due to out of phase action of his control. 
Phugoid oscillation must be highly considered in take-offs 
and ITandings. Short period is uSually out of pilot control, 
but is of importance in automatic control. 

The requirements that an aircraft should satisfy are 
written in appropriate manuals. It is usually imposed that 
the oscillation should not take more than one cycle to damp 


to half amplitude. 
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Flight test can be carried out to determine the 
oscillations characteristics under different conditions of 
Flight. Gyroscopic instruments can be placed and recordings 
of the longitudinal variables can be made possible to obtain. 
The aircraft is trimmed to the selected speed and height and 
longitudinal oscillation is initiated by an abrupt elevator 


deflection. 


4.3 PHUGOID AND SHORT PERIOD APPROXIMATIONS 
The fact that the angle of attack and hence, the 

w-velocity remains nearly unchanged in the phugoid mode, gave 
the idea of determining the phugoid quadratic by: 

(a) eliminating the w terms. 

(b) eliminating the pitching moment equation since no 
considerable pitching moment variation occurs, and 

(c) letting q=s9, obtaining thus the determinant: 


s- Xv 4 
-Zv -Uos 








(4-43) 


The characteristic phugoid quadratic comes out to be: 


2 d 2 
S Sarre =O (4-44) 


0 


For example, given, the characteristic quadratic is 
Ss +6.0097s + 0.09047 = B (4-45) 
which has roots 


-8.6049 + j 9.9681 (4-46) 





These roots are approximately equal to the exact phugoid 
FOOtES. 

From Equation 4-44 it is clear that the natural frequency 
of the phugoid is heavily dependent upon the 2 derivative 
while the damping of the phugoid is dependent upon the vu 
derivative. It is also obvious that high speeds tend to 
reduce the phugoid frequency and to increase the 
corresponding period. 

Since for subsonic straight level flight the 2 
derivative is simply expressed as: 

Zy = -29/ Vo (4-47) 


the following approximate relations hold for the phugoid 


Wry = 12 =. (4-48) 
a fe. Do Xv 
pH” oy q (4-49) 


Aircraft designers, in order to achieve larger values of 
phugoid damping, must increase the value of the *v derivative 
but this derivative is proportional to the drag. 
Fortunately, in take-offs and landings where phugoid motion 
1s undesired, high drag devices such as flaps are used 
anyway. 

A more precise approximation for the phugoid quadratic is 
given by leaving the w terms and the Mu term which in modern 
crafts is never zero. This is a three degree of freedom 


phugoid approximation. 
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s- Xv Xw q 
-Zv s-Zw Pus (4~-50) 
- Mv -My 0 


The characteristic phugoid quadratic becomes then 
1 i a g M 
2 v 

S +| Av 25 ae s- = (1. - qa 26] (4-51) 


which gives 


S My Ie 
Wine * | = (Lo - ae Zw) | (4-52) 


and 
Mu (Uo Xw -g) 


Ue Mw 


g My /2 (4-53) 
af = (2o- aie 2n)| 


Short period approximations can be made possible by: 
(a) eliminating the v-equations Since very small 
v-velocity change occurs. 
(b) setting all X-force derivatives to zero since 
very small X-force variation occurs. 
(c) letting q=s9, obtaining thus the determinant: 
-s-Zw -UVo 
~sMw-Mw s-Mq (4—-54) 
The characteristic short period quadratic comes out to 
be: 


52 + (-UM -Mq -Zw)s + (-U,Mw +My Zw) =O (4-55) 


G7 








For the example given, the characteristic quadratic is: 
s* + 4,.208s + 18.2556 = @ (4-56) 

which has roots 
=—2elimamt 7} 93.7187 (4-57) 

These roots are approximately equal to the exact short 
period roots. 

In summary, for typical flight conditions, the phugoid 
mode can be considered to consist of changes in the 
U-velocity and the 9-pitch angle, while the short period mode 
Can be considered to consist of changes in the w-velocity 
(angle of attack) and the pitch rate q. 

For the short period mode, the following approximate 
relations hold: 


Al2 


Wns, = (MgaZw -UoMw (4-58) 


-(U, Mw +Mq + Zw} 
ad 2 (Ma Zw ~ VoMw)"? (4-59) 
A critical condition occurs when the center of gravity is 


in that position for which 


Mg Zw = Uo Mw (4-69) 
For this case the characteristic equation becomes: 
s| s+ (-UMw -Ma -Zy)] =0 (4-61) 
which has roots 
0) Se UoMw + Mq + Zw (4-62) 


The non-zero root is usually negative which yields in a 


non-oscillatory and heavily damped exponential response. 
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4.4 SOLUTION TO THE LONGITUDINAL EQUATIONS 

The longitudinal equations form a set of linear first 
Order differential equations with constant coefficients. 
They are accompanied by a set of equal numbered initial 
conditions, 1.e. initial deviations from the equilibrium 
Flight conditions, and/or a number of control inputs. 

Many ways have been developed for solving sets of 
differential equations. Laplace Transform method is usually 
the simplest although it is not suitable for computer 
solutions. 

In the development of the solution, we would like to 
rearrange the equations in the matrix form: 

x (4) = A x (t) + Bult) (4-63) 
where x(t) is a n-dimensional time varying column vector 
named state vector, u(t) is a r-dimensional time varying 
column vector named control input vector, A 1s a nxn 
constant matrix and B is a nxr constant matrix. 


Taking the Laplace Transform of Equation 4-63 according 


to theorem 3 of Appendix A, we obtain: 


sX(s) = x) + AXG6) +B UG) (4-64) 
or 
(sI -A) X() = x) + BUG) (4-65) 
Or 
X(s) = (s1-A) x(0) + (sI-A) 8 U(s) (4-66) 


where IJ is the n x n identity matrix and x(@8) is the 


initial condition column vector. 
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The first term of the transformed solution corresponds to 
the initial conditions and is SSWABSGEen cots solution and 
the second term corresponds to the control input and is the 
particular solution. 

By inverse Laplace Transforming Equation 4-66 we obtain 
the time solution: 

x(t) =[2°Cs1-A)° | x0) + £4 [(sI-AY' Buc] (4-67) 

To implement this solution, in a digital computer, time 
domain formulation is followed. 

The key is to multiply the matrix differential 


Equation 4-63 by exp (-At) obtaining thus: 


ee _A 
e | x(t)-Ax@]-¢ te ut) Gis 
Or 
d , -At _ At 
xl? x 4) ae Put (4-69) 
Integrating between 9 and t, Equation 4-69 gives 
At t Ap 
e” xt) = x00) + | e Buct)dt ‘ieee 
Or 
At t  ACt-1) 
x (t) =@~ x(0) a g~ B u(t)dt (4-71) 


0 


The exponential term is expanded in Taylor series as: 

At i y 
er =[T+At Po Le bee AY (4-72) 
In reality, a series truncated to ten or twenty terms 1s 


often quite sufficient for normal accuracy requirements. 
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The computer program determines the values of the state 
variables at time intervals selected by the user. 
We now go back to the untransformed longitudinal 


equations and try to arrange them in the form of the matrix 


Equation 4-63. 
For convenience, 


replicated below: 


the untransformed equations are 


b - %)v - Ow) w 4 g@ = (X98 (4-73) 
AZo)v + w -(Za)w -(U)q = (258 (4-74) 
-(Mo) ~(Mw)w - (Mw) w + q -(Ma)q = (Ns )o (4-75) 
a ial (4-76) 
Equation 4-73 is written: 
POG Oe) wv - (g)@ BOs \S (Hag) 
Equation 4-74 is written: 
W = (2v)u +(Zw)w +CUo)q +(25)6 (4-78) 
Equation 4-75 is written: 
q = (Mv)v +(Mw)w +CMa)w +(Mq)q +WU5)5 (4-79) 
or substituting Equation 4-78 for w we obtain: 
a =(MytMiZv)v +CMw + Mi Zw) w + (Mio) +(Mg+ MZ) (4-88) 
Or 1£ we define: 
My = Mo + MwZo (4-81) 
Mw = Mwt Mw Zw (4-82) 
Mg = Mg + Mw Us (4-83) 
Ms* = M5 + Uw Zs (4-84) 
we will have: 
4 = (MB)u #CMB)w +(M$)q + (M5) 3 (4-85 ) 


le: 





It is further desirable to include the height of the 
. vehicle as another variable of interest. The equation for 
the height is obtained from the original Z-equation 
Equation 2-134 which for X% =@ becomes 
—7 = w -Vog (4-86) 
The left-hand side of this equation is just the z-axis 
acceleration, i.e. 
-d; Be hele w- Uo q (4-87) 
Integrating the last equation with respect to time we 
obtain 
h = -w+0,0 (4-88) 
the constant of integration being the initial w-velocity 
perturbation which is zero. Equation 4-88 describes the 
height perturbation around the nominal value. 
Collecting Equations 4-77, 4-78, 4-85 and 4-88 in matrix 
form we obtain the so called state variable representation of 
the longitudinal dynamics model, the state variables being, 


U,WrQ, 0, and h. 


v Xv Gr 0 -4 oO v XS 

W zy Zw Up O O W Z5 

q |= M*y Mt, MG 0 O q | + M‘S| 9 

2) O O { O O 2) O (4-89) 
h O -4 0 Use, On} | O 


ee 





State variable representation form is perfectly suitable 
for computer solution. Many ordinary differential equation 
solver software routines are available for solution. Most of 
these are provided with plotting capabilities of the state 
variables. 

The solution that follows corresponds to the given 
example with initial conditions and control inputs indicated 
in the summary of specifications tables. 

It can be seen from the plots that the v-velocity, the 
9-pitch angle and the h-height are heavily excited in the 
Phugoid motion, while the w-velocity and the q-pitch rate are 
heavily excited in the short period motion. 


For the phugoid approximation the following model may be 


used: 
v Xv “3 OQ U X§ 
© {=| -Zu/U. O Sener lezs/pal a (4-98) 
a O Vo O a O 


while for the short period approximation: 


De U, 
Me Ma. 


Ww w 25 
q , 4 M*s 


4.5 LONGITUDINAL AERODYNAMIC TRANSFER FUNCTIONS 


i d (4-91) 























Up to now we went through the analysis of the dynamical 
behavior of the flight vehicle, by making use of the 
mathematical model that adequately characterizes the 


longitudinal characteristic motions. In any dynamical system 


ly 3 





two types of variables exist, the input and the output 
variables. The input variables influence the output 
variables in a manner determined by the dynamics of the 
system. 

The objective of this section is that of determining the 
nature of those input =- output relations called transfer 
functions. 

For the longitudinal model, the input variable was the 
elevator control deflection $4 and perhaps the initial 
conditions or disturbances that can be treated as inputs, 
while the output variables were v,w,q,9 and h. 

The corresponding transfer functions are obtained by 
solving simultaneously at a time, the transformed equations 
of motion for the output variable of interest, keeping all 
initial conditions zero. 


The transfer function 9(s)/o0(s) is 


s- Xu - Xw O x5 

-Zy S- 2w U5 

-Mw -sMw-Mw s-Mq M9 
Xs) Z . PE 'ty fo s?+ Bos + Co 
WW) Aste Bake Cots DsvE Ashe Ble Gg Deve (4797) 

where 

Ag = Ms +MwZs (4-93) 
Bo = Xs5(Z,Miw +Mv) + Z5 (Mw -XuMw) - M3 (Xv+Zw) (4-94) 
Cg = X5CZy Mw -ZwMu)+ 25 (MuXw-MwXv) + M5 (Zw Xv -XwZv) = (4-95) 
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Or 


OG) Ag (5+ Z0,) Cs+ Zee) 
SG) Asts Bs2t Cs2 +Ds +E (4-96) 


where %,, 2, are the roots of the numerator polynomial 
Called zeroes of the transfer function. 


Similarily we can derive 


48) _ s6@) Ag s (s+Ze,) (s+ Ze.) 





As) &) Ast, Bo3+Cs2e DstE (a=ne 
ws) Aw s?+ Bw $?24Gws + Dw 
ats) Acts Bs3+Cs2+Ds+ & (4-98) 
where 
ON ee AS (4-99) 
Bw = X§lv -Z3 (Xo + Mg) + USU (4-198) 
Cw = X83 CU, Mp “Zu Mq ) + 28 Xv Mg -UoM Xv (4-161) 
Dw = q ( 23 My - M3SZv) (4-192) 


This transfer function may take the following forms 





wes) Aw (st¢2wi) (57+ 25wWw5S 4 Ww ) 

o(s) Ac4 4 Bs? 4 CsteDs+E (4-193) 
Or 

w(s) Aw (st2w,.) CS+2w,) Cs+Zws) 

&@) Ast 4 Bs?+Cs2eDs +E (4-104) 


We can therefore have three zeroes, two complex and one 
real or three real. 


Other longitudinal transfer functions are: 


Xs ) An ae For > GES ads 
as) Ag 4 Bg? + Cot*+Ds+ E (4-195) 


ILS 





where 
ae = X§ 
Bo = -X5 C Zw + Mg +UoMw) t Z5 Xw 


Ga XS C Zw Mg - Vo Mw) - 25 (Xw Mg taMw) + U5 (UoXw-9) 
Dp = g CMs Zw -25 Mw) 


This transfer function may also take the forms 


v5) Tey (CR chy, (PS SS aia 








dG) 7 As4 L BsF4 Co2+Ds+ [= 
Or 
vls) Av (St Zon) (s+Zv, ) Gcz Ze 
Ss) Acts Bs3+Cs2+Ds +E 
Similarily 
h(s) Hse as + Chs + Dp, 
§(s) s ( Ast + BsCs22Ds +E) 
where 


A, =~ 25 

Bh = - X8Zu +25 (Mg +loMw + Xv) 

Ch = X85 2v CMqg+lb Mw) -~25[¥o( Ma ¢Uo Mir) - UoMy] - USVo Mw 

Dn = - Ms Uo( 2w Mv -MwZv) + Z5 [ Mu (UoXw-g) -U, Mw Xo] 
+ M5 [Ur Zw Xv - Zv CUo%w -g) ] 


and may take for forms 


he) An Cst%,)(s°+ 25,95 + wp ) 
us) S(Acd + Be2+Cs2+ Ds+E) 


Or 


h(s) An (s4Zp,) Cst+Zy, )(s+2h3) 
ats) s (Ast + Bs? + Go? + Ds+E ) 


vas 


(4-196) 
(4-107) 
(4-198) 
(4-199) 


(4-119) 


(4-111) 


(4-112) 


(4-113) 


(4-114) 


(4-115) 


(4-116) 


(4-117) 


(4-118) 





Another transfer function of interest is that relating 


the z-axis acceleration to the control input, i.e. 


als)  -sth@) An s CStZp,) C 57+ 25,5 +H% ) 





ws) 0té‘ésGS)! Ae Dans (4-119) 
Or 
QAz(s) a Ah =) (st2p,) Cst Zh ) Cs+2p; ) 
A(s) As4+ Bs? +Cs24 Ds +E 
(4-129) 


Those transfer functions are used in automatic or manual 
Flight control systems, Since they specify the input - output 
relationships of the variables. For example, to find the 
open loop response 9(t) due to a unit step elevator input we 


use the transfer function Equation 4-96 


06) :| =| Ge te s (s+Z, )(s+Z5, ) 56) 


as) Ast 4 Bet+Cs24 Dstt | 


(4-121) 
Taking the inverse Laplace Transform of Equation 4-121 we 


obtain the time response. 


4.6 LONGITUDINAL EQUATIONS IN NON-DIMENSIONAL SYSTEMS 

The equations of motion of a flight vehicle, can be 
Written in a number of different forms, depending on the 
particular axes system selected for the definition of 
variables and the stability analysis. 

We have just described the body axes system in the 
dimensional form, in which most stability work is carried out 


fer. 1). 


a7 





In most aerodynamic problems it iS more convenient to use 
non-dimensional coefficients to represent the parameters 
involved. By this non-dimensional system, the equations of 
motion can be described in the body axes or wind axes system 
as well. 

The non-dimensional coefficients have the advantage that 
the effect of speed, size and air density are automatically 
accounted. However, it is required to express the mass and 
the time in non-dimensional forms as well. Therefore, one 
cannot obtain real time solutions. 

In the non-dimensional body axes system the longitudinal 


equations of motion are given as reference [Ref 2). 


dps - Cx, ~Cy, Cie v O 
2C,, - Cz, aus-s Cz. = Ge 7 dus => Cg “| & f= G:. aA8 (4-122) 
-Cm, ~§ Cm, -Cmy Lg Go SoM 0 S (m; + Cmny 


where the non-dimensional mass, inertia and time are 
expressed as: 

(A = 7m /9 Se (4-123) 
2 being the characteristic length taken as the half of the 


aerodynamic chord 


i.e MY 
: p se? (4-124) 
area (4-125) 
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In the non-dimensional wind axes system the longitudinal 


equations are [(Ref. 3] 


Sst G. 5 (@y-G) 5% 1) O 
CL s+ Cle =Q male | can m (4-126) 
0 sCm, + Cmy sCmg-s'Iy 110 | 1sCm,+Cma 


where the non-dimensional mass, inertia and time are 


expressed as: 





p= m/pSée (4-127) 
+ Syy /m 
-—Yyi = =? 

pc (4-128) 

c= am /'p Sus (4-129) 

4% t/z (4-13) 


Was; 





TABLE 4-1 


SPECIFICATIONS FOR COMPUTER SOLUTION OF THE 
LONGITUDINAL EQUATIONS DEDICATED FOR PHUGOID 
MOTION DUE TO INITIAL CONDITIONS 


Variables and Initial Conditions 

u = 5.9 ft/sec 

w= 2.5 ft/sec 

q = 86.95 rad/sec 

8 = 9.075 rad 

h=190.90 ft 

t = 0.9 

Stability Derivatives and Constants 

Xv = —-9.8097 Zv = -@.9955 Mu 

Xw = 9.0016 Zw = -1.43 Mw 

Xs = 0.0 Z5 = ~69.8 Mw 
Mg 
Ms 

Special Functions De 


My = Mv +My Zo 
a = Mw + My Lw 


Ma = Ma + Uw Vo 


Derivatives 

v= Xyv +Xww -gO + Xd 

A PAU me SPAY tVoq + 258 

q = Mov+ Maw+Mag +S d 
6 

h 


=- 


~ 4 
= -W+ Vo 


Outputs 


- g = 32.174 
-8.9235 U = 668.9 
-8.9913 
-1.92 
=26 .10 


Control Input 
a= 9.8 


U, Wr 9,9, hvs time at interval 9.9625 


end calculation when time > 2598.9 
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Figure 4-5 


U-Velocity vs Time for Specifications of Table 4-1 
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Figure 4-7 


gq-Pitch Rate vs Time for Specifications of Table 4-l 
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h-Height vs Time for Specifications of Table 4-l 
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TABLE 4-2 


SPECIFICATIONS FOR COMPUTER SOLUTION OF THE 
LONGITUDINAL EQUATIONS DEDICATED FOR SHORT PERIOD 
MOTION DUE TO INITIAL CONDITIONS 


6.0 qp= Ss2eln4 
-0.02355 v,= 660.9 
-9 .0013 

-1.92 

~26 .19 


Variables and Initial Conditions 

vu = 5.9 

w= 2.5 

gq = 0.05 

8 = 8.975 

h = 10.0 

t = 9.0 

Stability Derivatives and Constants 

X = —-@.0097 Zu = -0.0955 Mu 

Xw = 0.0016 Zw = -1.43 Mw 

Xj = 6.0 Zs = ~69.8 Ma 
He 
Ms 

Special Functions De 


Mo = MyeMwZy 


icw= Mews Mb-z,, 


+ , 

. = Mq + Ma Vo 
ME = M5 + MwZs 
Derivatives 


[wath 
( 


Xp v + Xww -g8 t X50 


w= Zypv + Zww +Uog + 255 

q = Mi y+ Miw + Mag + 39 
fee Hire 6 

6 = 4 

Outputs 


Control Input 
>= 0.0 


Ur, Wr Ge9, hvs time at interval 0.80125 


end calculation when time 75.90 
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Figure 4-19 
u-Velocity vs Time for Specifications of Table 4-2 
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Figure 4-11 


w-Velocity vs Time for Specifications of Table 4-2 
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Figure 4-12 


q-Pitch Rate vs Time for Specifications of Table 4-2 
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Figure 4-13 


Q@-Pitch Angle vs Time for Specifications of Table 4-2 
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Figure 4-14 


h-Height vs Time for Specifications of Table 4-2 
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TABLE 4-3 


SPECIFICATIONS FOR COMPUTER SOLUTION OF THE 
LONGITUDINAL EQUATIONS DUE TO INITIAL 
CONDITIONS AND CONTROL INPUTS 


Variables and Initial Conditions 

UV = 3. 

w= 2.5 

gq = 09.05 

8 = 9.075 

h = 10.90 

t = 9.0 

Stability Derivatives and Constants 

Xv = -9.0097 Zo = -8@.9955 Mo = @. g = 32.174 
Xw = 9.0016 Zw = —-1.43 My = -@.9235 Y. = 660.9 
X5 = 9.0 Zi; = -69.8 My, = -0.0013 

Ma = -1.92 
Special Functions D. Control Input 
d= 9.01 for 8<t<10.90 

. 6= -8.01 for 10.9<t<20.8 
Mv = Wy + Uw Zo 6= 9.0 for t>20.0 

Mw = Mw +Uy Zy 

+ : 

Mq = Mg + Mi Uo 

MZ = U5 + Uw 23 

Derivatives 
v= Xv v + Xww -g9 + X$0 
w= Zuv + Zww + Log + 259 
4 - MZu+ Wow + Ma 4 + MS 8 
ht = w+ U0 
Qe 4 
Outputs 


U, Wr, G,90, h vs time at interval 9.0625 


end calculation when time > 259.9 
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v-Velocity vs Time for Specifications of Table 4-3 
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Figure 4-16 


w-Velocity vs Time for Specifications of Table 4-3 
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Figure 4-17 


q-Pitch Rate vs Time for Specifications of Table 4-3 
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Figure 4-18 


©-Pitch Angle vs Time for Specifications of Table 4-3 
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h-Height vs Time for Specifications of Table 4-3 
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CHAPTER 5 


uh E AM 


5.1 INTRODUCTION 

In this chapter we will analyze the lateral directional 
motion of a conventional airplance, by folliwng the same 
treatment as in the preceding chapter on longitudinal 
dynamics. 

In the lateral directional case also, any particular 
response will be associated as in the longitudinal case, by 
its corresponding characteristic root location in the complex 
plane. 

Further simplified set of equations will be developed 
that apply to the individual modes of motion. Finally, 
computer solutions will be emphasized as well. 

As an example, the lateral-directional motion of the same 
aA€lrcraft as in the preceding chapter will be analyzed with 


the lateral-directional parameters indicated below: 


iOS 





altitude (ft) 20.000 
weight (lb) 30.500 
mach mumber 0.638 
air speed (f{thec) G60 


Yy = - 0.0829 lg = -4.77 Ng = 3.55 
Ys = 7.656 Lo = - 4.695 Np = -0,0025 
ly = OAT76 Nr = -0.0957 
le = 7.25 Ne = -06I5 
le = 066 Ne = -4.383 
Ixz/yy = 0.0663 Txz/t, = 0.0370 


5.2 CHARACTERISTIC POLYNOMIAL OF LATERAL—DIRECTIONAL 

MOTION 

Lateral-directional equations of motion are further 
Simplified by neglecting the effect of Yy, Yp, Yr, Ly, and Ny 
derivatives which are of no importance in lateral directional 
dynamics. Furthermore, a straight level flight is assumed as 
the reference flight condition (i.e. y=08). Hence, 


Equation 2-153 is reduced to: 


3-Yy ) { 3p o|| Bl | Ye 
-le S- -S es -Ly 0) O p LS 
- Ne ~$ FE -Np s-Np O O rl=| Ns{ 9 ai 

O ! O -5 b O 

: 0 O t 0 -§ yp O 


ILS) 











The characteristic polynomial, as known, is found by 
evaluation of the above determinant which comes out to be: 
As* + ee Cee a De +E =O (eZ) 
where 
Ie 


A 7 Ixx Jzz (5-3) 





ee -Y, (1- 2s) -1, - ne - Sy - 221, 





Ixx Izz Ixx izz (5-4) 
C= Ne + bp Y#Nr) + Np (22 W-Le} + Yo(e le Ne} BZ lp gy 
T, 
D = -Nglp + Yy ( NpLy -LoNr) t NpLg - Z (i+ ra N@ Ng) (5-6) 
g 
pos (Le Nr - Nal) (5-7) 


Equation 5-2 is referred to as the lateral-directional 
“~eapplity quartic. The roots in this case for any 
conventional vehicle are two real negatives or positives 
which are associated with two convergent or divergent motions 
and one pair of complex roots associated with an oscillatory 
motion. 

The overall lateral directional motion of the vehicle is 
a Superposition or a combination of all three modes. 

The complex pair of roots correspond to the as named 
Dutch Roll oscillation, while the large real root corresponds 
to the Roll Subsidence mode and the small root to the Spiral 


mode. 


200 





For the numerical example given, we find: 


A = 9.9975 _ (5-8) 
B = 1.867 (5-9) 
C = 3.6841 (5-10) 
D = 6.2637 (5-11) 
E = -0.0086 (5-12) 


and the characteristic polynomial is: 
0.9975 s* + 1.86753 4 3.G84I s? + 6.2657 5 - 00085 =O (Sade: ) 


which has roots: 


-9.0465 + 31.8784 (5-14) 
+0.014 (5-15) 
-1.7861 (5-16) 


The pair of roots, ~9.9465 + j1.8784, correspond to the 
Dutch Roll mode, while the root +@.8@014 corresponds to the 
Spiral mode and the root -1.7801 to the Roll Subsidence mode. 

Further, it is more convenient to put the characteristic 
polynomial in the form: 

(s7+25,, wry, $ +Wn_, )(s- 1; )(S-res ) (5-17) 
so that the natural frequency, the damping and the real roots 
are obtained by inspection while the damped frequency of the 
Dutch Roll can be determined by: 

Wye = Woe (1- She) (518) 

The root plotting in the complex plane will appear as in 


Figure 5-l. 


Z0 1. 








Figure 5-l 


Lateral-Directional Roots and Characteristic 
Quantities in the Complex Plane 


The imaginary roots of the characteristic polynomial can 
be written: 
“Spq Wien tj Wing (Bag) 
Or 
‘ 2 
~ Sop Wing + { Wry. ({- Soe) (5-26) 
Other quantities of importance in the lateral-directional 
dynamics are: the Dutch Roll period denoted by Tye given by: 
Noe = 2r/uy (5-21) 
the half amplitude for the Dutch Roll given by: 


0.6934 


ta 5p ~ Sa Wing (52995 
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the half amplitude for the Roll Subsidence mode (if stable) 


0,693 
Ves (5-23) 
and the half amplitude for the Spiral mode (if stable) 


ty, Rs 


0.6954 


yn oe 
M2 rs US OY.) 


For our example, the characteristic polynomial of 


Equation 5-13 can be factored in the form: 


(5? + 0.0935 + 3.5305) ( s- 0.0014) (s+ 4.7801) (5-25) 
From this by inspection we can determine: 

Wop = ( 3.5305)" « 1879 ~~ rad/sec (5-26) 

Sop = 0,093 /2- 1819 = 0.0247 (5-27) 

ae = 0.6931 /0.0247- 4.879 = 14.9051 sec (5-28) 

4ig,. = 0.6931 /;yg0; = 0.3894 sec (5-29) 


Many of the basic ideas involving longitudinal stability, 
also apply to lateral and directional stability. Lateral and 
directional stability are interrelated due to cross coupling 
effects of roll and yaw motion, that give rise to the three 
modes of motion discussed so far. 

In the Spiral mode, the airframe seems directionally 
nearly stable, with no sideslip but it is excited in a banked 
turn. The unstable mode is found usually in large finned 
airplanes with no dihedral, so that the side forces developed 
tend to turn the plane into the relative wind but since the 
outer wing travels faster, generates more lift and the plane 


will roll to still a higher bank angle. The flight of the 
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airplane is a heavy banked turn of ever-decreasing radius, 
i.e. a tightening spiral. The stable mode is found in 
airplanes with dihedral where the flight is a slight banked 
turn of gradually increasing radius. The spiral mode is 


illustrated in Figure 5-2. 


spiral convergence 


spiral divergence 


Figure 5-2 


Spiral Stable and Unstable Modes 


In the Roll Subsistence mode the airframe seems laterally 
nearly stable with no bank but it is excited in a sideslip 
turn. The unstable mode is present when side forces generate 
yawing moments that continue to increase the initial sideslip 
disturbance. This condition may continue until the airplane 
is broadside to the relative wind. Airplanes with successive 
weathercock stability, usually have this mode stable where an 
initial sideslip tends to die out. The Roll Subsistence mode 


is illustrated in Figure 5-3. 
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Figure 5-3 


Roll Subsistence Stable and Unstable Modes 


Dutch Roll is a motion which has the characteristics of 
the Spiral convergent and Roll Subsistence divergent mode. 
It is found in airplanes with strong lateral stability but 
weak directional. If a sideslip disturbance occurs and the 
plane is yawing in one direction, rolling moments are 
developed in a counter motion, i.e. yaw to the left is 
CauSing roll to the right and the plane wags its tail from 
Side to side. Dutch Roll mode is illustrated in Figure 5-4. 

The damped frequency of the Dutch Roll is indicated by 
the imaginary part of the complex root, while the real part 
indicates the damping. To reduce the effects of Dutch Roll, 
an increase in the directional and a decrease in the lateral 


Stability is required. 
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Figure 5-4 


Dutch Roll Mode 


5.3 DUTCH ROLL, SPIRAL, AND ROLL SUBSISTENCE MODE 

APPROXIMATIONS 

By analogy with the longitudinal short period and phugoid 
approximations, we can reduce the degrees of freedom to 
Obtain Dutch Roll Spiral and Roll Subsistence mode 
approximations. 

Dutch Roll is a motion which is characterized by 
negligible variations of roll and sum of the rolling moments 
nearly zero. The Dutch Roll quadratic is obtained from the 
matrix Equation 5-1, by: (a) eliminating the 9, and p 
equations, and (b) setting all rolling moment derivatives to 
zero obtaining thus, the determinant: 

5- Yv 4 
-Ng s-Nr (5-39) 
The characteristic Dutch Roll quadratic comes out to be: 


st? -(YW+tNr)s + (YvNr+Ng) =0 (5-31) 
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For the example given, the characteristic Dutch Roll 

quadratic is 

57 40.1786s + 3.5579 =0 (5-32) 
which has the roots 

-9.0893 + j1.8841 (5-33) 

These roots approximate the exact Dutch Roll roots. 

From Equation 5-31 it is clear that the frequency of the 
Dutch Roll is heavily dependent on the weathercock derivative 
Né and the damping on the yaw damping derivative Nr and the 
Yy derivative. 

All the derivatives involved in the Dutch Roll mode are 
dependent mainly from the vertical tail, and secondarily, on 
the wings and fuselage. Thus, proper care has to be taken 
For the vertical tail as far as the Dutch Roll 
Specifications. It is usually required that Dutch Roll 
oscillations should be damped to half amplitude within one 
cycle. 

A first approximation to the Spiral and Roll Subsistence 
mode is obtained by expanding the determinant, Equation 5-39, 
by the roll damping term (s-Lp) or s(s-Ljp. The determinant 


takes the form of: 


ey O { 
O s(s-Lp) O (5-34) 
Na 0 s-Ny 
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The characteristic equation comes out to be: 
s(s-Lp) | s?- (W+Nr)s + ( WNr+Ne)] =O (5-35) 

Here, the free s corresponds to the Spiral mode, the 
(s-Lp) corresponds to the Roll Subsistence mode, and the 
quadratic to the Dutch Roll mode. 

By this approximation, for the given example, the Spiral 
root is zero, the Roll Subsistence root is -1.695 and the 
Dutch Roll mode as previously -@.9893 + j1.8841. 

A more precise approximation for the Spiral and Roll 
Subsistence modes, is obtained by expanding the determinant 

O - We { 
—le s(s-Lp)  -br (5-36) 
-Ne s-Np -S-Nr 


This lead to 


Ne s°+(NpLe - Nelp - Le) 8 : : (LeNn- Nebr) =0 (5~37) 
Or 
t+ (Het 1,-22)s. 5 (4-1) =0 (5-38) 
For the given example one obtain the characteristic 
polynomial 
s*+1.7639 s -0.0024 =0 (5-39) 
which has roots 
@.0014, -1.7653 (5-49) 


These roots approximate the exact Spiral and Roll 


Subsistence roots. 
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5.4 SOLUTION TO THE LATERAL<DIRECTIONAL EQUATIONS 
Following the same approach as in the preceeding chapter, 
we first try to arrange the lateral~directional equations in 
the form 
x(t) = A x(t) +B uf) (5-41) 
For convenience, the untransformed equations are 


replicated below: 


B - (Yw)B +r -(3/b)o = (V¥s')8 (5-42) 
“le)B + 6 -Clp)p - (Tx2/tex) r - CLr)r = (18) 5 (5-43) 
~(Ne)B - (Tez/Tze) 6 - (Np) p +f - (Nr)r =(Na)S (5-44) 

6 =p (5-45) 

OS hy (5-46) 


Equation 5-42 is written 
B= (W)B -r +(4a)e + (Ys )5 Gan) 
Equation 5-43 is written 
) (La) B + Clp) p +(CIxzf,,)r +(Cly)r+(le)S +Us)5 (5-48) 
Equation 5-44 is written 
r = (Ne) a (Ixz/t,, ) ? + (Ne) p t(Ne)r +CNe)& HN5)S (5-49) 
In Equations 5-47, 5-48 and 5-49, the control inputs are 
expanded in terms of aileron (§) and rudder (3%) deflections, 
as indicated by Equation 3-143 neglecting the dotted 


derivatives lze,4e . 
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Substituting Equation 5-49 into Equation 5-48 we obtain 


pet Yn Me Let Ya Neg | Le if Ne 
Taf, Stet | Pe/tonte 


ls + dee Ixx Ne r Ne ee Nc 


4 - 1? x2 /TyyT22 4- echt 
(5-50) 


Similarly substituting Equation 5-48 into Equation 5-49 


we obtain 


ye Ne thef Lp ‘ Np + Ae Lp i Nr + we lr f 
\ - Lc en ee 4 - Iie les Ae ee ieedibes 


2 Ne phe. le . Ne plah, ls P 
4. 17) a Vis T*xz AxxTzz 
(5-51) 


We now define for convenience 


a hat aoe Né 


} = a ee 
5-52 
F I- Vee Atos 


ibs - Ly tte Np 


ae 5-53 
‘hs I ‘ee/ TyxIzz 


etc., so that Equation 5-58 and 5-51 can be written 
p= (le) 2 +(e) (er +e SE +E) 5 (5-54) 
r= (NE) B +(Np)p + (NE) r + WEE + (NE )S a 


PAY 





Collecting Equations 5-47, 5-54, 5-55, 5-45 and 5-46 in 
matrix form we obtain the state variable representation of 
the lateral directional dynamics model, the state variables 


being 4, pr r, ¢ andy. 


‘oe We -| sh OF oo || B Caeea(c 
6 ep ot © oo Ojlo' Ik & : 
r |= Ng No Nv O O O rj+| Ne Ne | | (5-56) 
d O 4 fo) 0 O O ® O 0 
| yp | 0 { O O O We : ef Bo : 


The solution that follows corresponds to the given 
example with initial conditions and control inputs indicated 
in the summary of specifications tables. 


For the Dutch Roll approximation the following model may 


be used 
B Yv -{ O p O Ye 
rj= Ne Nr OM Cole Ne Ne | : | (5-57) 
y O 4 O y O O 


which includes all the parameters of the determinant, 


Equation 5-30, while for the Roll Subsistence and Spiral 


approximation 
3 oO 8060 yi! jh oO || 6 O Ye 
Pale Lp Le O O}| p de ls : 
gee P= Ne No Ny O OF liye alee Ne Nx | | (5-58) 
p fe) { O O O fi) fo) O 
w O 0 O O wp | | O 0 | 


Ze 





which includes all the parameters of the determinant, 


Equation 5-36. 


5.5 LATERAL-DIRECTIONAL AERODYNAMIC TRANSFER FUNCTIONS 
Following a same approach as in the preceeding chapter, 
we can develope lateral-directional aerodynamix transfer 
functions. 
For the lateral-directional model, the input variables 
were the aileron control deflection (§) and the rudder 


control deflection (5) while the output variables were (@ p, ©, 


¢andv. 
Thus 
Pg 
Ys 0 ! -9/v, ) 
Ls s-Lp -s YY =r oO O 
Ns -s 3 -Np s-Nr O fe) 
0 4 0 ~$S Oo 
0 





My 


36) 5 | (Noe + Bs7+ Cs2iDs +t] 


_ Ag s?+Bas?+Gs+ Dz 


- As++Bs*+ Cs2+Ds +E (5-59) 


where 
Ag a ve Gia Vy2/TyeIzz ) (5-60) 
Be = -YS [Lp Ne eG Np eG, Le] - B34 Ls Hs (5~61) 
C6 is Y's CLp Nr - Np Ly) t Ng lp ~L5Np th, (lg +t Ne) (5-62) 
Da = Yh, (NsLr -LgNr) (5-63) 
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This transfer function usually takes the form 


Bo) _ Ag (s+Zg,)(s+Z, ) (s+Zg,) 
S56) As4.Bs?+Cs2+Ds+€ (5-64) 


Other lateral-directional transfer functions of interest 


are 
As) hy 5? + Bos + Cy 
Ss) Ag44 Boke Cs24Ds +E (5-65) 
where 
hy = 5 + ty NE Re) 
By = Ys (Le + ™ Ay Ne)-Ls (Net) © Ns (Le- 5% Y) (5-67) 
Co = YE CleNe Ug Nr) + La(WNeeNp)- Ns Cle + Yobe) (5-68) 


and takes the form 


ds) Ab (52+ 254095 +P ) 





Ie Ea ee (5-69) 
Also 
rs) Ar s?+ Brs*+Gs +Dr 
5s) Asts Bs3+Cs*+ Dstt (5-76) 
where 
Ar = NS+# tA Ls (5=71) 
Be = YS (Ne + PeA,, Lg) +L5 (Np-7*74,. Yv) - Ns (Yy +Lp) (5-72) 
Cr = YS (UeNp -Nglp) -L5YWwNp + Ng YL (5-73) 
Dr = Wy, CLsNg - Nols) | Sy) 
and takes the following forms 
rs) Ay Cst25) (32+ 25, s + wr ) 
als) Ast. Bs3+ Cs2eDs+€ (5-75) 


ZS 





Or 


ris) Ae (S#2r,) (s+2e,) (s#2r,) 
$0) Aste Bs?+ Cs24 DS t+ E (5-76) 





It is further desirable to include the transfer function 
relating the acceleration on the y-axis to the control input. 
The acceleration on the y-axis is given from the original 
V-equation, Equation 2-135, which for =§ becomes 

dy = Vv t Vor - 9b (SI 


and the cor responding transfer function is: 





Ay (S) _ Ag s4 + Ba s+, s°+Das + Ea 

aS) As4+ Bs? + Cs*4 Ds¢E (5-78) 
where 
Aa = Ys (A- Tx /ty Too ) (5-79) 
Ba = =Y5 lp te + Np + Bef, Lp) (5-80) 
Ca = Y8( No +lpNr -Nplr + 5/42, Le) ~Voy (74, Ls + Ns) (5-81) 
Da =Ys[Npla-Z le -Ne (lp +%, care) + Vou] Ls (2 - Np) +No ~ i +p) ee 
Ea = a [Ye (lee -Nelr) + Yv (NS Le -L5 Nr)| (5-83) 
and takes the form 

Ay(s) Aa ($+Za,) (s+2a, )(s+2a,) Cs+Za,) 

36) Ast4 Bs'+ Co24Ds +E (5-84) 


5.6 LATERAL-DIRECTIONAL EQUATIONS IN NON-DIMENSIONAL 
SYSTEMS 
In the non-dimensional body axes system the 


lateral-directional equations of motion are [Ref. 2]: 
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dys -Cy, -sQy,-G, Qu -Cy, Bi|o & 


5 
- Cog $7 La -SCp, -Ste -Co. |-] @ |= Q, @, iq (oz 65)) 
S) 
- Cng cot. -SOny Sle aiCye r Cn Cn 


where the non-dimensional mass, inertia terms and time are 


expressed as 


y= m/oSe (5-86) 
@ being now the half of the wing span (f= b/2 ) 
by a leo (5-87) 
Le = Iez/pSe’ (5-88) 
Le = Ixz/ose% (BSS) 
47 = l/y, (5-90) 


In the non-dimensional wind axes system the lateral 


directional equations are [Ref. 3]: 


-25 +Cy, Zs C, B 6 Gy, 
Ho, OC 3s@ -Skx |.) yisi@ || (5-91) 
p Gn, = Cn, -sl, ; s(Mp> 4) Cre Cie 


where the non-dimensional mass, inertia terms and time are 


expressed as 


u =m loSb (5-92) 
e/a 

Ix = 2 ne (5-93) 
lz2/m 

12 = ub? (5-94) 

a m /0 SV. (5-95) 

ey (5-96) 
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TABLE 5-1 


SPECIFICATIONS FOR COMPUTER SOLUTION OF THE LATERAL- 


DIRECTIONAL EQUATIONS DUE TO INITIAL CONDITIONS 


A. Variables and Initial Conditions 
6 = 6.85 
P = 6.05 
r = 6.075 
b= 6.1 
V= 6.2 
t = 0.0 
B. Stability Derviatives and Constants 
Yy = -9@.0829 Lg = -4.77 Ne = 3.55 Ixz/Ixx = 0.0663 
Ys = 7.656 Lp = -1.695 Np = -@.0025 Ixz/Izz = 6.03786 
Lr = 6.1776 Nr = -@.9957 g = 32.714 
Le = 27.25 Ny = -@.615 Y = 660.0 
Ls = 9.666 Ns = -1.383 
C. Special Functions 
ie Js SOD 
Oo = ee (elec avy ee | 
a 
la = L lp +0 2/tyy) Ne] /D Ng = [Ng + (he ) Le] /D 
Lp = | Lp + ( Txz/ q, x) No | /D > = | Np +(Ixz/q22 ) Lo] /D 
ag = | Ly + CLee/Tyx) Nr] /D NY > | Ny +( dxe/7,2) Lr J /D 
We = [ de + (te2/z,,) Ne] /D E= [Ne +(e/tz2) Le /D 
ie = ) Le ~ i277, | Ns ]/2 z = | Ng xz) ) hes 7D. 
D. Control Inputs 
5= 8.6 
>= 6.0 
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TABLE 5-1 (Continued) 


E. Derivatives 


)) 


Vy B -r +(9/U0)b + Ye 5 
le ( + lbp ro a a +le +i 
Ne + Nop + Nrr + NEE +Ng3 


€. 2. a oO WD 
y 


F. Outputs 
, Pr Ce >, Y VS time at interval 9.0125 


end calculation when time > 59.@ 
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Figure 5-5 


&-Sideslip Angle vs Time for Specifications of Table 5-1 


218 





————— 
——— + 


— 


13 

| 

ei 

ie Oram 


| 
ag 
| 


| 
Ie 








ia — Cra aa Gh 





2452 


Figure 5-6 
p~Roll Rate vs Time for Specifications of Table 5-l 
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r~Yaw~Rate vs Time for Specifications of Table 5-1 
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Figure 5-8 


$-Bank Angle vs Time for Specifications of Table 5-l 
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Figure 5-9 


w-Yaw Angle vs Time for Specifications of Table 5-l 
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TABLE 5-2 


SPECIFICATIONS FOR COMPUTER SOLUTION OF THE LATERAL- 
DIRECTIONAL EQUATIONS DUE TO INITIAL CONDITIONS AND 
CONTROL INPUTS 


Ng 


= 855 Ixz/Ixx = 
= -§.0025 Ixz/Izz = 
= -§ .8957 es 32.07 4 
= -9.615 Vo = 660.90 
= -1.383 

Na =LNe + (Txz/7,, la] ID 
Ny = [ Np t(Ixz/tz2 )lpJ /? 
Né = [ Nr + CIxz/Izz) Lr] /D 


NE = [ Ne +(Ixz/Iez2)Lgj /D 
Ne = [ Ng +(Tx2/322) Lg J/D 


A. Variables and Initial Conditions 
8 = 9.95 
p = @.95 
r = 9.9075 
6 = 8.1 
w= 8.2 
t = 0.0 
Be. Stability Derivatives and Constants 
Yv = -8@.0829 Le = -4.77 
Lr = 9.1776 
Meee) 29 425 
Ls = 0.666 
C. Special Functions 
Ve = ‘x [Uo 
a Tenia: ) 
Ue = [le + ( Iv2/J,,) Ne | /D 
Up ~ | bp + Clut/t,x ) Np | ID 
be = | ly ¢ ( Ix2/4,x) Nr J /D 
le = Chey + (1x2/q..) Ng] /D 
le = | Lz + (Lee JIxx) Ng jJ/D 
D. Control Inputs 


Uy Ut UTS 


Om Ol for @<t<2.5 
-9.01 for 2.5<t<5.90 
Umue ton to} .8 
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vw ut vt UT 
uoeu we ul 
Qi ana 
e Qe e 


foe) 0<t< 15.8 


0 
itor wlo.0 <tc lie 5 


® for 520.0 


9.0663 
9.0370 


ieeOGe 7 ost <2 0 
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TABLE 5-2 (Continued) 


Derivatives 

A= WR -r + Ole)o +55 d=9 
pe eto eltr s+ Ee els b = 
r - NEB tNSp eNEr e NEE + Ns 

Outputs 


, Pr Yr >, vs time at interval @.6125 
end calculation when time > 50.0 
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p-Roll Rate vs Time for Specifications of Table 5-2 
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Figure 5-10 


8-Sideslip Angle vs Time for Specifications of Table 5-2 
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Figure 5-12 
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r-Yaw Rate vs Time for Specifications of Table 5 
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Figure 5-13 


od-Bank Angle vs Time for Specifications of Table 5-2 
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Figure 5-14 


WY-Yaw Angle vs Time for Specifications of Table 5-2 
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CHAPTER 6 


STATE VARIABLE FEEDBACK LAW 


6.1 INTRODUCTION 

The automatic-flight control functions of an aircraft are 
provided by an automatic flight control system which is 
referred to as autopilot. Autopilot statilizes the dynamic 
response of the aircraft providing the best response to 
disturbances and control inputs. Today, almost all types of 
alrcrafts, helicopters and missiles are equipped with 
autopilots with complication depending on the particular 
application. 

The equipment comprising the autopilot includes sensors, 
controllers, estimators, computer and actuators. To damp the 
vehicle responses to disturbances and control inputs, the 
rates of the angular motion in roll, pitch and yaw are sensed 
by gyroscopic instruments and fed back to the flight control 
computer. The computer combines the feedback information 
with the desired pilot response and computes commands which 
are applied to the actuators that deflect the control 
Surfaces. 

In this chapter we will discuss the state variable 


feedback method as a means to control the vehicle. 
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6.2 THE FEEDBACK CONCEPT 

miicemMecatrve feedback attempts to modify thie 
characteristic equation of the system so that one can get the 
Specified by the requirements response. 

Consider as an example that the roll motion of an 
aircraft is described by the equation 
© = lop + lsd (6-1) 


@ 


The open loop transfer function is 


PS) Ls 
as) S- Lp (6-2) 


Suppose that the p-variable is continuously measured and 


fed back as shown in Figure 6-l. 


SR 
SS A A 





Fygure 651 


Simple Feedback System for Roll Stabilization 


Zell 





The closed loop transfer function is 





PO as /Cs-bp) Ls 
os) A+ Kp Li/cs-lp) S$ + (hls -lp) (6-3) 


It is clear that the single characteristic root changed 
from Lo to Kplgy-Lp. Usually Lp is negative and Ls positive. 
So the system will have a faster response indicated by the 
greater negative root of the characteristic polynomial. 

In automatic control systems all the variables of the 
system are fed back with corresponding gains K's chosen to 
give a particular desired response. 

In general, the system is described in matrix form as 

e=Ax+By Slee 
with initial condition vector 
X ©) = Xo (6-5) 
Applying the feedback law 
ue-K x +g Nee?) 
where 
Ko= PRR, K, <.. 9%) (o— 7) 
is the feedback gain matrix and ug the desired input, one 


obtains the system 


x = (A- 


(@ 


Keay 6 ua (6-8) 
The matrix A-BK now determines the new characteristic 
roots of the system. Note that the characteristic roots of 


the system, Equation 6-4, are obtained by the determinant 


[sI -Al =O (6-9) 


Ze 





while for the system, Equation 6-8, by the determinant 
|s1-(A-8K)I <0 (6-10) 

Since the response of the system is determined by the 
location of the characteristic roots, any desired response 
will specify the characteristic polynomial of the system 
which if equated with Equation 6-19 the feedback gains are 
obtained. 

Feeding back all the variables means that somehow all of 
them must be measured or estimated. Usually not all of the 
variables are meaSured but the ones that are not measured are 
estimated. For example, the pitch angle can be obtained by 
integration of the q-pitch rate. 

O == 7 (6-11) 

Even unstable systems can be stabilized with state 
variable feedback provided that the system is controlable. 

Most times the feedback gains are selected such that for 
the system to meet a performance criterion. In this case 
optical control methods are used to select the gains. 

In the next pages a computer solution is shown of the 
given example for the longitudinal and lateral directional 
case uSing the state variable feedback concept. 
Specifications for computer solution are shown in Tables 6-l 


and 6-2. 
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TABLE 6-1 


SPECIFICATIONS FOR COMPUTER SOLUTION OF THE LONGITUDINAL 
EQUATIONS DUE TO INITIAL CONDITIONS, CONTROL INPUTS AND 
STATE VARIABLE FEEDBACK STABILIZATION 


0.0 g = 32.174 
-9.0235 uv. = 660.0 
-9.6013 

=I G22 

-26 .10 


rol Inputs 

@.01+ £ for 0<t<10.0 
-9.O01+ £ for 10.9<t<20.0 
06.0 + £ for t>20.0 


Variables and Initial Conditions 

v = 5.0 ft/sec 

w = 2.5 ft/sec 

q = @.05 rad/sec 

© = §@.075 rad 

h=16.0 ft 

t = 0.0 

Stability Deriviatives and Constants 

Xy = -@.0097 Z = -@.0955 My = 

Aw = 9.0016 Zw = -1 .43 Mw = 

Xs = 6.0 Z, = -69.8 My = 
Ma = 
Ms = 

Special Functions D. Cont 

My = My +MwZu im 

Mw = Uw + Ui Zw f= 


Ua c Ma + Uw Vo 
Us’ = 5 + Mw Z5 


Derivatives 

v= Xvv + Xww - 90 + X53 
w= 2vu + Zww tiba + 23d 
q= US + Mw + Uh g + HES 
oS 
h 


= -W + U.O 


Outputs 


-0.005 uy -90.0005 w + 0.022 4 
+ 0.0035 9 + 0,.00000014 h 


U, We Gr 9, Hh vs time of interval 9.0625 


end calculation when time > 250.90 
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Figure 6-2 


U=-Velocity vs Time for Specifications of Table 6-l 
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Figure 6-3 


w-Velocity vs Time for Specifications of Table 6-1 


236 





. 
. 
. 
. 
. 
. * 
ne ee ee 


O10 


a ee 
. : = . 
° : 
. 
. 
. 
° 
ne 
. 
ry 
‘ 
r 
rf 
° 
me aa tt a a Oe 


ao a eee et em ee 
. 





ns en 
. 


YW VS TIME 


eo enee om 2 eEtOL UNITS INCH. 
ener Of-Oe UNITS 


INCH. 


Figure 6-4 


q-Pitch Rate vs Time for Specifications of Table 6-l 
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Figure 6-5 


©-Pitch Angle vs Time for Specifications of Table 6-1 
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h-Height vs Time for Specifications of Table 6-1] 
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TABLE 6-2 


SPECIFICATIONS FOR COMPUTER SOLUTION OF THE LATERAL- 
DIRECTIONAL EQUATIONS DUE TO INITIAL CONDITIONS, CONTROL 
INPUTS AND STATE VARIABLE FEEDBACK STABILIZATION 


A. Variables and Initial Conditions 
8 = @.05 rad 
p = 8.05 rad/sec 
r = 9.075 rad/sec 
od = 9.1 rad 
Y= 9.2 rad 
t = 8.0 sec 
B. Stability Derivatives and Constants 
Yy = -@.0829 Le = -4.77 Ne = 3.55 Ixz/Ixx = 0.0663 
Y, = 7.656 Lp = -1.695 Np = -0.0025 Ixz/Izz = 0.0370 
Lr = 9.1776 Ne = -9.0957 q = 32.174 
Le = 27.25 Ny = -@.615 Vv, = 660.9 
Ls = 9.666 N, = -1.383 
C. Special Functions 
Ys = Mee 05 
Ds i-( Ve / Ixx Tzz) 
-*K = 
Ue = [ Lg + Ctxz/T,4) Ne J /D Ne aN ke (Txz/t2. Le] /D 
C= [bp +Cbe/tyx) Np] /D No = [ Np + (Txz/te2) Lp 1 /D 
Be = [ br + Clee /qy,) Nr] /D NF = ( Nr +(Tx2/z,7) br J /D 
ly = i Le cect x) Ne] /D Ne = [ Ne ncz/ Teo be jf 
= aes pexz | ty. ) Ne jf Ng = | Weenie Joy) he 1) D 


D. Control Inputs 


Pr — Omolet ese tor 9<t<2.5 mee Oe or £ror, O0<t< 15.0 

§ = -@.01 + 5£ for 2.5<t<5.9 Peewee st For 15.,.0<t<1/.5 

E = 9.0 + Sf for t>5.9 S= -0.1 + SE for 17.5<t<20.0 
S$= 9.0 + S£ for t>20.0 

Ff = ~O0USBR -0.222 p -0,209r -0 324 @ -O2G67 pH 

sf = 


0.0515 B + 0.05069 + 0.4954 + 0.0932 +0169 H 
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E. 


TABLE 6-2 (Continued) 


Derivatives 

B=YWRetr + (dd +YZ8 

p= brbo rte + ees 
Pos Ng 6 +*Nop +Nrr s+ NEEB+NES 
¢ =? 

ber 

Outputs 


Br Pe Le by YW VS time at interval 6.0125 
end calculation when time > 5@.90 
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Figure 6-7 


@-Sideslip Angle vs Time for Specifications of Table 6-2 
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Figure 6-8 


p-Roll Rate vs Time for Specifications of Table 6-2 
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R VS TIME 
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Figure 6-9 


r-Yaw Rate vs Time for Specifications of Table 6-2 
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Figure 6-19 


?-Bank Angle vs Time for Specifications of Table 6-2 
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~-Yaw Angle vs Time for Specifications of Table 6-2 
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CHAPTER 7 


AT 0 S E 


7.1 $INTRODUCTION 

The prceeding chapters have been oriented primarily 
towards the dynamics of the conventional piloted airplanes. 
The equations of motion for the missile are still valid since 
it can be thought as an automatically controlled aircraft. 

Missiles are found in different geometry considerations. 
Cruise missiles usually have no rotational symmetry, i.e. 
rotation through 360° is required along the missile axis to 
leave the external shape unchanged. Configurations with 186° 
rotational symmetry are seldom encountered while 128° or 
higher-order rotational symmetry are more typical. 

In this chapter we will consider the cruciform missile, 
i.e., with 96° rotational symmetry to take the advantage of 
symmetries. This configuration is shown in Figure 7-l. This 
misSile may be considered as typical for today's large 
existing missiles. 

In the next sections we will consider the simplification 
to the equations of motion due to symmetry and the 
complication of concerning the coupling effect of 
longitudinal and lateral-directional motions. The last 
happens because strong rolling moments are produced when the 


missile is simultaneously pitched and yawed. 
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Figure 7-1 


Cruciform Missile Configuration 


For this particular model we will assume that control 
inputs are coming from the vertical and horizontal fin 
deflectors of the tail. A deflection of the vertical fins 
Stands for a rudder deflection while a deflection of the 
horizontal fins stands for an elevator deflection . These 
are the only existing control inputs. Roll control is 
obtained by differential deflection of one vertical or 


horizontal fin. 
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7.2 SYMMETRY CONSIDERATIONS 

For the missile's equations of motion, the reference body 
axes system shown in Figure 7-1 will be used as in the 
airplane case, and straight level horizontal flight will be 
assumed as the nominal flight condition. 

Considering the Semi symmecries: we find first that 
in addition to the aircraft symmetries the xy plane is a 
Plane of symmetry and that makes the inertia tensor of the 
missile having only diagonal terms and as a consequence the 
body axes coinsides with the principal axes. Further 
Iyy = I2zz. 

The fact that the missile is not designed as an airplane, 
gives no tendency to remain in the same roll orientation. 
Vertical and horizontal find deflections can be optimized to 
give zero bank angle in turning maneuvers. This condition is 
highly desirable in missiles as well as in airplanes because 
the lifting forces are most efficiently generated and 
lift-—to-drag ratio is maximun. 

However, accidential actuator error signals, asymmetrical 
loading of the control and lifting surfaces in Supersonic 
flight and atmospheric disturbances may introduce large roll 
rates and large roll rates cause cross-coupling between the 
longitudinal and the lateral-directional motion. This makes 
the equations of motion not separable in a longitudinal anda 
lateral-directional set which was found so convenient in the 


aircraft case. 
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We therefore find the equations of motion, Equation 2-134 
and 2-135, combined in one set of equations and modified as 


follows: 


4 ° 
mans =o +g (j=8 
{ 
Fry = W eG z (7-2) 
A M = } 
Tyy $ ane 
A : 
wri es V + Vor - 9 (7=h) 
{ _ 
ms (7-5) 
{ 
ic ee (7-6) 
with auxilliary relations 
p= ob (7-7) 
= (7-8) 
re w (7-9) 


Before expanding the forces and moments in terms of 
Stability derivatives we note tht due to symmetry the Y-force 
is the same as the Z=-force which makes v=w and 

Yv = Zw (7-18) 


Yo = Zw (7-11) 
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Also the N-yawing moment is the same as the M-pitching moment 


which makes r=-q and 


Yr =-Zq (7a) 
Nv = -Mw (7-13) 
Ny = -Mw (7-14) 
Nr = Mq (7-15) 


because Iyy = Izz 


In addition, for the control derivatives one has 


Ys = -Zn (7-16) 
Nz = Mn (7-17) 
Ni = Ma (7-18) 


7.3 EQUATIONS OF MOTION FOR THE NON-ROLLING MISSILE 
If we expand Equations 7-1 to 7-6 in terms of stability 


derivatives we obtain 


= Xvv + Xww -g9 + Xnn (7-19) 
Ww = 2v vp + Zw W + 224) + Zn Y) (7-28) 
q = Muv + Myw + Uy Ww +Mq q +Unn (7-21) 


Or substituting Equation 7-28 


q = (My+ Mw2v)yu + (My +MwZw)w + (Mq + Mw Uo) 4 +(Mn+Mw2n)n (7-22) 


me iy =a0er + ao + ¥ze5 (7-22) 
) = LyVv + Lo p tirr Fy oats (7-24) 
eos Nyy + Nvv + Npp + Nrr + NzS (7-25) 


At Equation 7-25 we' would like to keep the term N,v 
because it is the same in magnitude as M.* Substituting 


Equation 7-23 for v 


Zonk 





r= (NvtNvY )v +Npp +(Nr-NiUo)r +gNid + (Ne+Nv Ye )5 (7-26) 
EE Ene missive is not rolling the equation is 
eliminated as well as the p-stability derivatives. Also, if 
the v-velocity perturbation is assumed zero, the equation 
is eliminated as well as the v~-Stability derivatives. 


Therefore, one obtains the following model for the 


non-rolling missile 


w-= Zy W + Uo g + Znn (7-27) 
G2 (Mw+UwZw) w + (Mas Mog +CMn tM Zn) n (7-28) 
y = WY -Vor + Y23 (7-29) 
r= CNvtNi Yv)v +(Nr-Nvo dr 4 (Ne + Ny )3 (7-30) 


Or more conveniently in matrix form 


w Zw IDI: O O W an @ 
; ¥ * + 
Oo O O 
5 Ae Lh Be : eon 
v O 0 Yv -U, Vv O Xz, S 








f atta e O Ny Ne dlyvl}] oo Ne 


This model may be augmented with the auxiliary equations 
and transfer functions and state variable feedback can be 


applied for proper stabilization. 
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Theorem l: 
Theorem 2: 
Theorem 3: 
Theorem 4: 
Theorem 5: 


Theorem 6: 


APPENDIX A 


SOME LAPLACE TRANSFORMS 


f(t) 4970 LS Fk)} = FG) 
9 (4) 4 
u(t) (unit step) 1s 
r(t) or t Cunit ramp ) 4/52 
+? 2/s3 
grat ‘/ (sta) 
os ‘/ C seat)? 
cos wt s/ (s?4w?) 
sin wt w/( s?+ar) 
2 E** cos Cwt +@) Ce { Desc 


SOME LAPLACE TRANSFORM THEOREMS 
Lfasebaayt = af) +b RG) 
LS «fait = « FG) 
Li df /dt } = s FG)- fe) 
Ld Fe) /c? } = Fey - sfe) - del, 
L{[ A L(t) /He} = s* FG) - sf) -s didn |, - df /at?|, 
£§ (° pide} = FOL 


iT 
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